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An  experimental  study  has  been  conducted  to  examine  the 
three-dimensional  natural  convection  heat  transfer  from  an 
array  of  simulated  electronic  components  immersed  in  a 
chamber  filled  with  Fluorinert  FC-75,  a  commercially 
available  dielectric  liquid.  The  top  and  bottom  walls  of 
the  chamber  were  maintained  at  uniform  temperature  while  all 
other  surfaces  were  insulated.  The  simulated  components 
were  in  the  form  of  a  3  x  3  array  of  discrete  protruding 
aluminum  blocks,  each  with  geometrical  dimensions  of  a  20 
Pin  Dual-inline-Package.  The  components  were  electrically 
powered  resulting  in  a  range  of  energy  dissipation  levels 
from  0.1  to  3.1  watts.  Flow  visualization  in  steady  state 
was  accomplished  using  Magnesium  particles  illuminated  by  a 
Helium  Neon  laser  plane.  Component  surface  temperature 
measurements  allowed  determination  of  the  heat  transfer 
characteristics.  Timewise  fluctuations  of  temperature  at 
several  locations  were  measured  with  increasing  power 
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I.  INTRODUCTION 

A.  STATEMENT  OF  PROBLEM 

As  microelectronic  components  continue  to  decrease  in 
size  and  increase  in  component  density,  the  problem  of  heat 
dissipation  becomes  increasingly  more  critical  in  the  design 
of  circuits  and  components.  A  f ive-order-of-magnitude 
increase  in  circuit  integration  in  the  last  25  years  has 
been  associated  with  successive  revolutions  in  components 
technology. 

With  the  recent  development  of  CMOS,  industry  has  seen 
an  almost  three-order-of-magnitude  reduction  in  size,  an 
order-of-magnitude  increase  in  the  characteristic  component 
dimensions  and  perhaps,  most  importantly,  a  precipitous  drop 
in  transition  switching  energy  from  10-9  J  in  1960  to  nearly 
10”13  J  in  present  day  devices.  It  is  anticipated  that 
power  dissipation  of  100  W  in  a  10  million  component  chip 
could  be  encountered  by  1990.  This  would  result  in 
component  heat  fluxes  in  excess  of  100  W/cm2  [Bar-Cohen, 
1987,  pp.  159-175].  The  strong  requirement  to  keep 
components  at  low  temperature  levels,  100°C  to  125°C,  does 
exist,  since  for  every  20°C  decrease  in  components 
temperature,  the  components  failure  rates  are  cut  in  half 
[Oktay,  1986,  pp.  36-42]. 
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B.  BACKGROUND 

Air  cooling  is  the  standard  design  in  today's  systems. 
Advances  in  air  cooling  technology  now  include  impingement 
cooling  and  combined  conduction  and  air  cooling  technology 
which  offers  even  higher  cooling  rates  than  natural  or 
forced  convection  air  cooling.  However,  air  as  a  heat 
transfer  medium  is  extremely  limited  not  only  in  its  ability 
to  transfer  heat  but  also  in  its  capacity  to  carry  the 
transferred  energy  with  a  reasonable  temperature  increase 
[Baker,  1972,  pp.  213-222]. 

Special  component  surfaces  have  been  investigated. 
Features  such  as  Turbulence  Promoting  Fins  (TPF)  have 
significantly  higher  cooling  enhancements  when  used  with 
special  gases  such  as  Helium  as  the  closed  cycle  coolant 
[Kishimoto  et  al.,  1983,  pp.  286-93].  By  doing  so,  a  300% 
greater  dissipation  rate  can  be  reached  than  in  an  open 
cycle  with  air  as  the  coolant. 

Forced  convection  in  liquids  offers  greater  heat 
transfer  coefficients  than  does  natural  convection  in 
liquids  or  air.  However,  forced  convection  introduces  two 
major  problems:  (i)  greater  amounts  of  expensive  coolant 
are  required  and  (ii)  pumps  and  external  plumbing  increase 
the  size  and  reduce  the  reliability  by  introducing  the 
possibility  of  vibration  related  failures  to  the  system 
[Baker,  1973,  pp.  163-173]. 


Boiling  can  produce  very  large  scale  heat  transfer 
coefficients,  but  it  also  introduces  major  problems:  (i) 
boiling  restricts  the  physical  design  of  the  system,  (ii)  if 
the  surface  is  very  smooth,  as  most  electronic  components 
are,  a  high  degree  of  superheat  may  be  required  so  that 
nucleation  can  occur.  The  occurrence  of  periodic  boiling 
due  to  the  high  local  heat  flux  may  result  in  severe  local 
turbulence  and  mechanical  damage  and  (iii)  boiling  may 
result  in  the  concentration  of  impurities  at  the  nucleation 
site  [Baker,  1973,  pp.  163-173]. 

C.  STUDIES  ON  NATURAL  CONVECTION  COOLING  OF  ELECTRONIC 

DEVICES 

1 .  Numerical  Studies 

Of  the  different  heat  transfer  modes  mentioned 
above,  immersion  cooling  by  natural  convection  offers  not 
only  high  heat  transfer  potential,  but  also  a  practical 
means  of  controlling  circuit  temperatures.  However,  only 
limited  heat  transfer  data  are  available  for  immersion 
cooling  of  discrete,  protruding  heat  sources  [Liu  et  al., 
1987,  pp.  309-30].  Early  studies  of  natural  convection  heat 
transfer  were  motivated  by  applications  such  as  energy 
conservation  and  solar  energy  collection.  The  pertinent 
length  scales  were  several  order  of  magnitude  higher  then 
those  applicable  in  cooling  of  electronic  equipment. 
However,  recent  studies  have  focused  on  computational  and 


experimental  studies  of  natural  convection  cooling  in 
electronic  components. 

[Han  et  al.,  1984,  pp.  21-28]  have  used  numerical 
simulation  to  investigate  the  transient  laminar  natural 
convection  in  a  rectangular  enclosure  with  a  localized 
constant  heat  flux  source.  [Acharya  et  al.,  1983,  pp.  855- 
66]  have  investigated  the  natural  convection  phenomenon  in 
an  inclined  square  box  shaped  enclosure  which  has  energy 
sources  in  it.  [Jaluria,  1985,  pp.  883-92]  investigated  the 
interacting  natural  convection  flows  generated  by  flush, 
isolated  thermal  energy  sources,  such  as  electronic 
components,  located  on  a  vertical  adiabatic  surface. 
[Jaluria,  1982,  pp.  223-27]  also  analytically  studied  the 
natural  convection  flow  due  to  multiple  isolated  heated 
elements  located  on  a  vertical  surface.  [Torok,  1984,  pp. 
49-57]  used  the  finite  element  technique  to  predict 
component  and  board  temperatures  and  heat  fluxes. 

[Kuhn  et  al.,  1986,  pp.  55-62]  have  investigated  the 
transient  three-dimensional  natural  convective  flow  and  heat 
transfer  within  a  rectangular  box  following  a  sudden 
temperature  change  on  a  rectangular  portion  of  one  vertical 
wall.  They  found  that  when  two  elements  are  symmetrically 
placed  on  the  wall,  the  distance  between  them  has  only  a 
slight  effect  on  the  mean  heat  transfer  distribution. 
However,  it  does  have  a  significant  effect  on  the  local 
Nusselt  number.  [Shakerin  et  al.,  1986,  pp.  17-22] 


investigated  the  natural  convection  flow  adjacent  to  a 
heated  wall  with  single  and  repeated,  two-dimensional, 
rectangular  roughness  elements.  [Oosthuizen  et  al.,  1987, 
pp.  1-7]  studied  the  heat  transfer  rates  from  an  element  of 
square  cross-sectional  shape  held  at  a  uniform  temperature 
and  mounted  on  a  vertical  adiabatic  wall  of  a  tilted  square 
enclosure.  They  found  that  the  mean  heat  transfer  rate  from 
the  block  did  not  vary  continuously  with  the  angle  of 
inclination. 

[Liu  et  al.,  1987,  pp.  486-97]  have  conducted  a 
finite  difference  numerical  study  of  natural  convection  flow 
in  a  rectangular  enclosure  of  an  array  of  simulated 
components  that  are  mounted  on  a  vertical  wall.  The  results 
of  this  numerical  study  have  been  presented  for  enclosure 
widths  of  18  and  30  mm.  It  was  concluded  that: 

-  The  temperature  field  in  the  enclosure  is  characterized 
by  boundary  layer  regions  surrounding  individual 
components . 

-  There  are  only  small  interferences  among  the  components, 
especially  for  the  lower  rows. 

-  The  maximum  temperatures  on  the  components  surfaces  are 
located  on  the  upper  horizontal  surfaces  of  the 
components.  Lower  temperatures  have  been  found  on  the 
lower  horizontal  faces  of  the  components  in  the  bottom 
row. 

-  The  results  for  the  30  mm  width  are  essentially  the  same 
as  those  for  the  18  mm,  indicating  only  a  minor 
dependence  on  the  enclosure  width. 

[Liu  et  al.,  1987,  pp.  309-30]  have  also  observed  a  local 
oscillatory  surface  temperature  response  in  the  same 


enclosure.  Their  study  found  that  the  components  surface 
temperatures  have  a  tendency  to  fluctuate  within  a  range  of 
approximately  3°C  with  a  period  of  4  seconds. 

2 .  Experimental  Studies 

Numerous  experimental  studies  have  been  conducted 
concurrently  with  the  numerical  studies  described  above. 
Two  of  the  earliest  studies  were  conducted  by  [Baker,  1972, 
pp.  213-22]  and  [Baker,  1973,  pp.  163-73].  The  purpose  of 
these  investigations  was  to  examine  liquid  immersion 
techniques  for  cooling  miniature  heat  sources,  specifically 
effects  of  component  size  on  the  convective  heat  transfer. 
[Goel  et  al.,  1986,  pp.  1341-46]  experimentally  investigated 
the  two-dimensional  natural  convection  flow  from  a  long, 
finite  size  heat  source  located  on  a  vertical  or  inclined 
plate  in  air.  The  principal  items  of  interest  were  the 
nonboundary  layer  effects  that  arise  in  the  flow, 
particularly  in  the  vicinity  of  the  source.  Flow  separation 
was  observed  at  large  inclinations. 

[Milanez  et  al.,  1986,  pp.  1347-52]  experimentally 
studied  the  interaction  of  two  horizontal,  wide  strips  on  a 
vertical  adiabatic  wall  in  both  air  and  water.  They  found 
the  results  of  their  experiment  in  good  agreement  with  past 
numerical  studies.  [Ortega  et  al.,  1986,  pp.  5-15] 
presented  the  results  for  a  series  of  experiments  in  air 
designed  to  measure  the  natural  convective  heat  transfer 
rates  and  the  thermal  characteristics  of  a  rectangular  array 


of  cubical  roughness  elements  mounted  on  an  insulated  plane 
wall,  with  and  without  a  shrouding  wall.  [Parks  et  al., 
1987,  pp.  91-109]  simulated  microelectronic  circuits  with 
thin  foil  heaters  for  both  flush  and  protruding  components. 
They  found  that  the  protruding  heaters  have  approximately 
15%  higher  coefficients  and  that  the  upper  heaters  have 
varying  coefficients  depending  upon  the  array  configuration. 

A  number  of  studies  have  been  conducted  at  the  Naval 
Postgraduate  school.  [Knock,  1983]  studied  the  effect  of 
the  location  of  a  single  protruding  heater  in  an  enclosure 
filled  with  water.  [Pamuk,  1987]  investigated  the  heat 
transfer  characteristics  of  an  array  of  nine  protruding 
components  immersed  in  a  dielectric  fluid  in  an  uninsulated 
chamber.  [Hazard,  1987]  studied  the  natural  convection 
liquid  cooling  of  a  vertically  arranged  components  array 
with  and  without  a  shrouding  wall. 

D.  OBJECTIVES 

This  experimental  work  was  motivated  by  the  numerical 
study  of  [Liu  et  al.,  1987,  pp.  486-97]  and  the  recent 
experimental  studies  of  [Pamuk,  1987].  The  objectives  for 
this  thesis  are: 

-  To  design  and  build  a  circuit  card  of  nine  simulated 
components  and  determine  the  rate  of  heat  transfer  for 
various  power  inputs. 

-  To  determine  the  natural  convection  flow  patterns 
developed  in  the  chamber  by  means  of  flow  visualization. 

-  To  correlate  the  data  in  terms  of  dimensionless  heat 
transfer  parameters  such  as  Rayleigh  vs.  Nusselt  number. 
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II.  EXPERIMENTAL  APPARATUS 


A.  GENERAL  CONSIDERATIONS 

This  experimental  study  is  a  continuation  of  past 
experiments  conducted  at  the  Naval  Postgraduate  School 
[Knock,  1983;  Pamuk,  1987;  Hazard,  1987].  The  purpose  of 
this  work  was  to  examine  natural  convection  heat  transfer 
and  fluid  flow  from  a  3  x  3  element  array  of  simulated 
electronic  components.  A  major  change  from  the  previous 
experiments  was  the  design  of  the  simulated  circuit  card. 
In  [Pamuk,  1987]  an  arrangement  of  vertical  components  was 
used.  The  present  study  utilized  a  horizontal  arrangement. 
A  decision  to  increase  the  card  thickness  was  made  to  reduce 
the  contamination  of  the  dielectric  fluid  by  the  loose 
insulating  material  previously  used  and  also  to  reduce  the 
conduction  loss  through  the  card. 

The  experimental  apparatus  consists  of  a  chamber  filled 
with  a  dielectric  fluid.  Top  view  of  the  chamber  is 
provided  in  Figure  2.1.  The  top  and  bottom  surfaces  of  the 
chamber  are  maintained  at  a  prescribed  temperature  by  means 
of  two  heat  exchangers  and  two  separate  chilled  water 
circulation  baths.  One  vertical  wall  of  the  chamber 
simulates  a  circuit  card.  It  contains  nine  simulated 
components  in  a  3  x  3  array,  while  all  other  vertical  walls 
are  insulated.  The  geometrical  arrangement  of  the  simulated 
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components  on  the  card  is  seen  in  Figure  2.2.  The 
components  were  chosen  to  simulate  an  array  of  20  pin  Dual¬ 
inline-packages  (DIPS) .  Their  dimensions  are  identical  to 
those  used  by  [Liu  et  al.,  1987,  pp.  486-97;  Liu  et  al.f 
1987,  pp.  309-30;  Pamuk,  1987]  in  order  to  permit  future 
comparisons  of  experiments  and  computations. 

The  heating  elements  were  made  of  aluminum  and  heated  by 
means  of  small  foil  type  resistance  heaters  attached  to  the 
bottom  surfaces.  The  temperatures  of  the  faces  of  the 
elements  were  measured  by  means  of  embedded  thermocouples. 
Data  were  recorded  and  processed  using  an  automatic  data 
acquisition  system.  The  major  components  of  the  system  are 
described  in  the  following  section. 


B.  COMPONENTS 

1.  Heating  Element 

Each  heating  element  was  a  block  of  aluminum  6mm  by 
8mm  by  24mm  (Figure  2.3).  The  temperatures  at  the  faces  of 
the  blocks  were  measured  by  means  of  thermocouples  placed 
within  grooves  (0.5mm  deep)  on  each  surface.  The  front 
surface  thermocouple  was  passed  through  the  block  by  means 
of  a  1.02mm  hole  drilled  through  the  block.  All  grooves 
were  filled  with  high  thermal  conductivity  epoxy  and 
smoothed  prior  to  attachment  to  the  simulated  circuit  card. 

2 .  Thermocouples 

All  thermocouples  were  10-mil  (0.254mm)  copper- 
constantan.  These  small  wires  were  chosen  to  minimize  the 
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material 

removed 

from  the 

block. 

The 

thermocouples 

were 

manufactured  by 

using  an 

arc  welder 

designed  for 

such 

purposes . 

The 

setting 

for  both 

current  time  and 

gas 

pressure  (argon)  was  kept  at  'medium. '  Due  to  the  extremely 
small  diameter  of  the  wire  the  gas  pressure  was  carefully 
regulated  to  less  than  5  psig.  The  thermocouples  were  also 
checked  for  continuity  due  to  the  fact  that  the  arc  welder 
tended  to  embrittle  the  copper  wire  at  the  bead. 

3 .  Heaters 

A  nearly  constant  heat  flux  condition  was  maintained 
at  the  base  of  the  block  by  use  of  thin  (0.18mm)  foil  type 
heaters  (Figure  2.4).  Each  of  the  heaters  was  measured  to 
have  a  resistance  between  10.5  and  10.6  ohms.  These  heaters 
were  designed  to  align  with  the  blocks.  Holes  were  provided 
within  the  heaters  to  pass  the  various  thermocouple  wires. 

The  power  leads  of  the  heaters  were  gold  plated  to 
ensure  a  low  electrical  contact  resistance.  Wires  from  the 
power  supply  were  attached  to  the  leads  by  means  of 
soldering.  After  ensuring  proper  alignment,  the  heaters 
were  bonded  to  the  blocks  by  means  of  a  high  thermal 
conductivity  epoxy. 

4 .  Simulated  Circuit  Card 

Three  grooves  of  0.5mm  depth  were  made  on  the  front 
of  the  card  in  order  to  allow  the  heaters  to  be  mounted 
below  the  surface  and  thus  prevent  the  dielectric  fluid  from 
seeping  under  the  blocks  and  introducing  a  variable,  unknown 
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resistance.  The  heating  elements  were  bonded  to  the  circuit 


card  with  a  low  thermal  conductivity  epoxy. 


5 .  Test  Chamber 


The  test  chamber  was  constructed  of  plexiglass  (13mm 


thick)  and  had  inside  dimensions  of  120mm  in  height,  144mm 


in  length,  and  30mm  in  width.  The  wires  were  fed  through 


the  grooves  in  the  back  of  the  circuit  card  to  a  "Tygon" 


tube  on  the  back  wall.  The  chamber  also  contained  two  small 


pressure  equalization  outlets,  located  at  the  top  and  the 


bottom  of  the  front  vertical  surface. 


Heat  Exchangers 


The  two  aluminum  heat  exchangers  measure  38mm  by 


65mm  by  274mm.  The  details  of  the  heat  exchangers  are 


available  in  [Knock,  1983],  Aluminum  plates  3mm  thick  were 


used  as  top  and  bottom  surfaces.  These  plates  were  grooved 


2.5mm  to  provide  slots  for  the  thermocouples.  This  left  a 


0.5mm  wall  thickness  of  aluminum  between  the  thermocouple 


bead  and  the  surface  in  contact  with  the  dielectric  fluid. 


The  resulting  conductive  thermal  resistance  was  found  to  be 


negligible.  The  heat  exchangers  were  each  attached  to 


individual  coolant  baths  by  means  of  "Tygon"  tubing.  Four 


10  mil  (0.254mm)  thermocouples,  two  on  the  top  and  two  on 


the  bottom,  were  used  to  read  the  heat  exchanger 


temperatures . 
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7 .  Flow  Visualization 


Visualization  was  accomplished  with  magnesium 
powder,  325  mesh.  Magnesium  has  a  density  of  1.92  g/cm3 
compared  to  1.7  6  g/cm3  for  FC-7  5,  resulting  in  a  large 
suspension  time  of  the  particles.  The  beams  from  the  light 
sources  were  spread  into  plane  sheets  which  illuminated  the 
particles  suspended  in  the  FC-75.  This  technique  allowed 
for  the  visualization  of  a  single  plane  of  the  flow  field. 
Other  planes  were  examined  by  the  realignment  of  the 
apparatus.  A  four  milli-watt  helium  neon  laser  and  a 
cylindrical  lens  were  used  for  the  time  lapse  camera  photos 
(Figure  2.5). 

8 .  Assembly 

After  sliding  the  card  into  the  slot  in  the  chamber, 
both  the  thermocouple  and  power  wires  were  run  out  through 
the  1.27cm  diameter  extension  pipe  on  the  back  of  the 
chamber  to  the  connected  "Tygon"  tubing  of  approximately  1.0 
meter  in  length. 

Rubber  o-rings  were  placed  in  grooves  located  on  the 
top  and  bottom  surfaces  of  the  chamber  and  the  thin  aluminum 
plates  were  set  in  place.  The  heat  exchangers  were  then  set 
on  the  plates  and  clamped  by  means  of  specially  made 
aluminum  clamps  to  ensure  a  proper  seal.  The  inlet  and 
outlet  tygon  tubes  were  attached  to  the  constant  temperature 
baths  for  each  heat  exchanger.  Finally,  the  pressure 
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equalization  tubes  were  attached  to  the  chamber, 
schematic  of  the  entire  assembly  is  shown  in  Figure  2.6. 
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III.  EXPERIMENTAL  PROCEDURE 

A.  APPARATUS  PREPARATION 

After  assembling  the  apparatus,  the  chamber  was  filled 
with  the  dielectric  fluid  through  the  top  vent  tube  by  means 
of  a  small  funnel  and  a  hypodermic  syringe.  To  damp 
disturbances,  the  entire  apparatus  was  placed  on  a  separate 
platform.  This  platform  also  served  as  a  means  of  leveling 
the  chamber  by  use  of  adjustable  screw  legs.  Since  the 
entire  apparatus  was  sealed  by  means  of  four  clamps  it  could 
be  tilted  in  order  to  remove  the  last  remaining  air  bubbles 
inside  the  chamber. 

The  use  of  a  single  constant  temperature  bath  providing 
cooling  water  in  series  to  the  top  and  bottom  heat 
exchangers  proved  inefficient  to  maintain  both  the  top  and 
the  bottom  heat  exchangers  at  the  same  constant  temperature. 
A  second  bath  was  added  so  that  each  heat  exchanger  was 
supplied  individually  and  could  be  adjusted  separately.  The 
ten  terminals  on  the  power  supply  panel  were  used  to  supply 
power  to  the  heaters  on  the  simulated  components.  The  first 
terminal  was  the  main  input,  the  next  nine  were  the 
terminals  to  the  individual  heaters  located  on  each 


component.  These  heaters  were  connected  in  series  with  the 
2  ohm  precision  resistors.  This  arrangement  was  repeated 


for  each  component.  These  nine  units  were  then  connected  in 
parallel  to  the  main  input. 

B.  DATA  ACQUISITION 

Prior  to  energizing  the  power  supply  or  the  constant 
temperature  baths  all  the  thermocouples  were  scanned  to 
measure  the  starting  steady  state  temperatures.  The 
constant  temperature  baths  were  set  at  10°C  by  means  of 
adjustment  dials  located  on  the  baths.  Throughout  each 
experiment  the  bath  supplying  the  top  heat  exchanger 
required  constant  monitoring.  As  the  power  levels 
increased,  the  bath  temperature  was  reduced  accordingly  to 
maintain  the  top  heat  exchanger  thermocouple  readings  at  a 
constant  10°C. 

The  desired  voltage  was  obtained  by  setting  the  power 
supply  to  the  desired  level.  The  reading  on  the  power 
supply  was  verified  by  means  of  reading  channel  number  61, 
the  main  power  input  channel,  on  the  data  acquisition  panel. 

Approximately  four  hours  were  allotted  for  the  system  to 
achieve  steady  state.  For  this  experiment,  the  steady  state 
condition  was  defined  as  a  majority  of  the  thermocouple 
readings  changing  less  than  0.2°C  during  three  separate 
readings  taken  at  ten  minute  intervals.  The  bath 
temperatures  were  checked  to  ensure  they  were  at  10°C.  If 
they  were  in  agreement  then  the  data  was  recorded  by  means 
of  the  program  "ACQUIRE”  (Appendix  C) .  This  program  scanned 
all  76  channels  and  converted  the  emf  signals  from  the 


thermocouples  to  temperatures  in  degrees  Celsius.  The  power 
to  the  heater  was  determined  by  utilizing  the  precision  2 
ohm  resistor.  This  resistor  was  connected  in  series  with 
the  heater  and  allowed  the  measurement  of  the  current 
through  this  configuration.  This  current  was  then 
multiplied  by  the  measured  voltage  of  the  heater  to  produce 
the  power  input  to  the  component. 

A  second  program  "Fastscan"  (Appendix  C)  was  developed 


to  scan  three  individual  channels.  Each  of  the  three 
channels,  numbers  1,  13,  32,  was  scanned  individually  at  a 
rate  of  3  scans  per  second  for  a  period  of  200  seconds.  The 
main  intention  in  doing  so  was  to  isolate  any  periodic 
behavior  caused  by  the  flow  past  these  three  thermocouples. 
A  third  program  "Plot”  (Appendix  C)  was  developed  to  plot 
these  results  on  a  HP  model  7470A  plotter. 

Flow  visualization  was  conducted  by  means  of  325  mesh 
Magnesium  particles  and  a  four  milliwatt  laser  focused  with 
a  cylindrical  lens.  The  magnesium  particles  were  injected 
into  the  chamber  by  way  of  the  bottom  vent  tube.  A  motor 
driven  Nikon  F3  camera  with  a  Micro-Nikkor  55mm  lens  and  a 


Nikon  MT-2  Intervalometer  was  used  for  the  long-time 
photography.  The  time  exposure  photographs  of  length  8  or 
20  seconds  were  taken  at  the  various  power  levels  where  heat 
transfer  data  had  been  obtained.  The  photographs  were  taken 
of  the  three  columns  of  components  by  placing  the  camera  at 
a  ninety  degree  angle  to  the  plane  created  by  the  laser. 


This  configuration  was  then  reversed  to  obtain  photographs 
of  three  separate  planes  in  the  chamber. 

C.  DATA  ANALYSIS 

1.  Determination  of  the  Nusselt  Number 

The  Nusselt  Number  is  defined  as  the  following: 

Nu  =  ip  (3.1) 

Kf 

where  L  is  a  characteristic  length,  taken  here  as  the  length 
of  the  component  in  the  vertical  direction,  kf  is  the 
thermal  conductivity  of  the  fluid,  and  h  is  the  convective 
heat  transfer  coefficient.  All  fluid  properties  were  taken 
to  be  a  function  of  the  average  film  temperature. 

To  evaluate  the  heat  loss,  one  thermocouple  was 
placed  directly  behind  each  of  the  heated  blocks.  A  series 
of  nine  thermocouples  were  also  attached  to  the  back  of  the 
card.  These  temperatures  were  used  to  calculate  the 
conduction  loss  through  the  back  of  the  card. 

One  dimensional  conduction  resistance,  Rcon,  for  the 
simulated  board  was  determined  from: 

Rcon  =  kA  (3.2) 

where  L  is  the  length  along  the  conduction  path  i.e.,  the 


thickness  of  the  board,  k  is  the  thermal  conductivity  of  the 
material  (Appendix  A)  and  A  is  the  area  normal  to  the 


conduction  path.  The  area  normal  to  the  path  of  conduction 
was  calculated  as  the  area  directly  behind  the  components 
back,  192mm2.  By  doing  so  spreading  losses  were  neglected. 
The  conduction  loss  is  then: 


T  ~  T, 
s  b 


(3.3) 


Power  input  to  each  component  is  calculated  from: 


<vin  -  VVh 


(3.4) 


where  V^n  is  the  input  voltage  applied  to  both  the  2  ohm 
resistor  and  the  heater  which  is  in  series  with  it.  Vjj  is 
the  voltage  drop  across  the  heater.  The  difference  between 
V^n  and  Vh  gives  the  voltage  drop  across  the  2  ohm  resistor. 

The  net  energy  added  to  the  fluid  in  the  chamber  is 
the  difference  between  Qin  and  Qcon  which  by  Newton's  Law  of 
cooling  is: 


Qconv  ~  Qnet  -  *lA(Tavg  “  Tc) 


“  Qin  “  Qcon 


(3.5) 


where , 


A  =  The  total  wetted  surface  area  of  the 
components , 

TaVg  =  Area  weighted  average  surface  temperature, 
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calculated  by, 


E  A .  T  . 

l  l 


(3.6) 


where , 


Temperature  of  the  i*-*1  surface  of  the 
components, 

Area  of  the  ith  surface  of  the  components. 
The  summation  of  the  index  i  varies  from 
1  to  5, 


and, 


Tc  =  Average  heat  exchanger  temperature. 


Thus, 


xconv _ 

Atot^Tavg  V 


(3.6) 


The  above  procedure  was  repeated  for  each  of  the 
nine  components  at  each  power  level.  The  Nusselt  Number  was 
then  calculated  from  Equation  3.1. 

2 .  Determination  of  the  Ravleiah  Numbers 

The  definition  of  the  temperature  based  Rayleigh 
number  is: 


Rat  =  Grt  *  Pr 


(3.7) 
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where  the  Grashof  Number  (Gr^)  is  defined  as 


and  the  Prandtl  Number  (Pr)  is  defined  as, 


Pr  =  - 

a 


(3.9) 


where , 


(pcp) 


All  the  properties  in  the  above  equations  are 
calculated  at  the  average  film  temperature: 


T  +  T 
Tf  =  av?-^  C 


(3.10) 


Therefore: 


Rat  = 


gSL  (T  -  T  ) 
^  ava  c 


(3.11) 


3 .  Determination  of  the  Flux  Based  Ravleiah  Number 

The  definition  c  *  the  flux  based  Rayleigh  number  is: 


Raf  =  Grf  *  Pr 


(3.12) 
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all  film  properties  are  calculated  at  the  average  film 
temperature  as  before  and: 


«i-40net 

Grf  =  2 

V  Atot 


(3.13) 


Functional  relationships  for  the  properties  at  the 
film  temperature  are  given  in  Appendix  A,  together  with 
sample  calculations.  Uncertainty  calculations  are  given  in 
Appendix  B. 
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IV.  RESULTS 


A.  FLOW  VISUALIZATION 

Flow  visualization  was  accomplished  using  Magnesium 
particles  as  a  flow  marker  and  a  Helium  Neon  laser  to 
illuminate  the  particles.  The  results  provided  insights 
into  the  flow  patterns  existing  at  various  power  dissipation 
levels.  As  expected,  the  particle  traces  suggest  highly 
three-dimensional  flow  patterns  at  all  power  levels. 
Particles  can  be  seen  entering  and  leaving  the  plane  of 
light  in  almost  all  the  long  time  exposure  photographs.  A 
multitude  of  flow  patterns  was  observed  for  the  wide  range 
of  input  power  levels  in  the  of  0.1  to  3.1  watts. 
Representative  flow  visualizations  are  presented  in  Figures 
4.1  to  4.14. 

1.  Flow  Patterns  at  Component  Midolanes 

The  following  flow  patterns  are  in  a  vertical  plane 
through  the  center  of  components  4,  5,  6  (Fig.  2.2).  The 
laser  was  placed  in  front  of  the  chamber  and  the  camera 
located  to  the  side  (Fig.  2.5).  The  camera  exposure  time 
was  set  at  10  seconds.  The  three  dark  regions  located  on 
the  left  side  of  the  photographs  are  the  three  rows  of 
components  as  seen  in  side  view.  The  top  and  bottom  of  the 
chamber  and  the  horizontal  walls  can  also  be  seen  in  the 
photos . 
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Figure  4.1  is  the  flow  pattern  for  a  power  input  of 


0.1  watts  and  indicates  a  pattern  that  contains  two 


individual  eddies  located  between  the  rows  of  components 


along  with  an  overall  circular  pattern  that  encompasses  the 


entire  chamber.  A  boundary  layer  like  flow  is  seen  in  the 


vicinity  of  the  heated  components.  A  dark  region  located  at 


the  bottom  of  the  photograph  is  evident  in  all  photos  taken 


throughout  the  experiment.  This  region  is  the  stagnant. 


dense  cold  fluid  layer  in  the  vicinity  of  the  bottom  heat 


exchanger.  This  region  proved  to  be  impermeable  at  all 


power  levels. 


Figure  4.2  for  power  input  of  0.3  watts  indicates 


the  same  overall  pattern  with  the  exception  of  a 


constriction  of  the  overall  circular  pattern  at  the 


approximate  height  of  the  second  row  of  components.  The 


slight  increase  in  the  upper  eddy  size  is  also  evident. 


Figure  4 . 3  for  a  power  input  of  1.5  watts  indicates 


an  increase  in  the  thickness  of  the  boundary  layer  emanating 


from  the  second  row  of  components  along  with  increased 


entrainment  of  the  particles  located  below  the  second  row  of 


components.  Also  evident  is  the  disappearance  of  the  two 


distinct  eddies  seen  earlier  and  strong  evidence  of  three 


dimensionality  in  the  flow  as  seen  by  the  crossing  of  the 


particle  paths. 
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Figure  4.4,  for  a  power  input  of  2.7  watts,  shows  a 
buoyant  fluid  layer  adjacent  to  the  components  and  random 


motion  throughout  the  remaining  chamber.  Also  seen  are 
entrainment  of  particles  at  the  approximate  height  of  the 
second  row  of  components,  and  evidence  of  an  area, 
immediately  above  the  black  region  at  the  photo's  bottom, 
where  particles  are  basically  suspended. 

Figure  4.5,  at  3.1  watts,  demonstrated  a  downward 
flow  pattern  in  the  vicinity  of  the  bottom  row.  This 
downward  flow  was  forced  past  the  face  of  the  bottom  chip  in 
this  column.  The  flow  also  demonstrated  a  perpendicular 
entrainment  of  particles  between  the  first  and  second  row  of 
components,  and  a  large  increase  in  the  boundary  layer  at 
the  third  row  of  components. 

Figures  4.6  and  4.7  show  flow  patterns  in  the 
vertical  plane  at  the  mid-span  of  components  7,  8  and  9. 

These  two  photographs  were  taken  approximately  5  seconds 
apart  at  a  power  setting  of  3.1  watts.  They  demonstrate  the 
strongly  time-dependent  behavior  of  the  flow.  In  Figure 
4.6,  entrainment  into  the  upward  moving  buoyant  layer  is 
seen  to  be  present  almost  as  far  down  as  the  bottom  of  the 
lowest  component.  Far  away  from  the  components  the  core 
fluid  exhibits  no  preferential  direction  of  motion.  In 
Figure  4.7  the  flow  pattern  has  changed.  Near  the  lowest 
row  of  components  the  flow  now  appears  almost  quiescent. 
Also,  the  buoyant  region  appears  thinner  than  in  Figure  4.6. 
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2 .  Flow  Patterns  in  Other  Planes 

Figure  4.8  is  a  schematic  of  the  three  planes  that 
were  investigated  when  the  positions  of  the  camera  and  the 
laser  were  reversed.  Figures  4.9,  4.10  and  4.11  are  at  a 
power  setting  of  0.1  watts  and  demonstrate  the  flow  patterns 
in  the  back,  middle  and  front  planes,  respectively.  In 


Figure  4.9,  the  back  plane,  the  flow  appears  to  be  somewhat 
like  a  large  circular  pattern  that  is  bounded  by  the 
physical  dimensions  of  the  chamber  and  exhibits  no 
noticeable  flow  directly  up  the  center  of  the  chamber. 
Figure  4.10,  the  middle  plane,  indicates  a  strong  upward 
flow  in  the  center  of  the  chamber  and  two  areas  of  eddy 
formation  located  approximately  mid-height  in  the  chamber. 
Finally,  Figure  4.11,  the  front  plane,  demonstrated 
approximately  the  same  behavior  as  the  middle  plane;  a 
general  area  of  motion  in  the  center  of  the  chamber  along 


Figure  4.8  Top  View  of  Planes 
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Figure  4.11  Front  Plane  @  0.1  Watts 
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Figures  4.12,  4.13  and  4.14  are  the  same  three 
planes  at  a  power  setting  of  3.1  watts.  Figure  4.12,  the 
back  plane,  indicates  regions  of  upward  flow  near  the 
components.  These  patterns  characterize  the  boundary  layer 
nature  of  transport  in  the  vicinity  of  the  elements. 
Relatively  short  particle  traces  in  regions  away  from  the 
components  appear  to  result  due  to  the  strongly  three- 
dimensional  transport.  The  dark  area  in  the  very  bottom  of 
the  photograph  is  again  pointed  out  to  demonstrate  the  cold 
isothermal  layer  developed  by  the  bottom  heat  exchanger. 

Figures  4 . 13  and  4 . 14  are  views  of  the  center  and 
front  planes  respectively.  The  flow  in  both  planes  appears 
strongly  time  dependent.  Figure  4.14  shows  relatively 
longer  particle  traces  than  Figure  4.13,  indicating  stronger 
three  dimensional  effects  near  the  components. 

These  13  photographs  illustrate  the  complex  three- 
dimensional,  time  dependent  flow  patterns  found  in  this 
study.  Obtaining  a  detailed  understanding  of  these  effects 
is  essential  in  achieving  predictive  capability  for  the  heat 
transfer  in  these  configurations. 

B.  REDUCTION  OF  HEAT  TRANSFER  DATA 

The  temperature  and  reduced  data  tables  corresponding  to 
the  eight  power  levels  are  as  follows:  (see  Figure  2.2  for 
identification  of  the  component  numbers) . 
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Figure  4.13  Middle  Plane  @  3.1  Watts 
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TABLE  1 

TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  0.3  W 


RESULTS  ARE  STORED  IN  FILE:  D80S8 !  T  30 


AMBIENT  TEMP  HAS :  23  C 

VOLTMETER  LEADIN',  -AS:  2.0332  *61 
BATH  T£MP  -AS:  t0  •' 

THIS  RON  HAS  RECORDED  HITH  FLOW  VIZ 


ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NOt : 

16.S4E+00 

!6.33E*00 

1 G . I SE  *00 

16. 1 6E  *00 

16. 1 1E+00 

1 6 . 96E  +  00 

POWER 

(  WA  T  TS ) : 

23.03E-02 

CHIP  N02: 

t  S .  B8E-00 

1 6 . 3  T  E  *00 

16.34E+00 

16. 19E*00 

16.30E+00 

1 7 . S3E+00 

POWER 

(WATTS) : 

29.84E-02 

CHIP  N03: 

16. 74E+00 

I6.32E*00 

I6.46E+00 

1 6 . 33E*00 

16.63E+Q0 

17.40E+00 

POWER 

(WATTS) : 

29.83E-02 

CHIP  N04: 

16.S2E+00 

16.23E*00 

1 6 . 08E+00 

IS . 82E*00 

1 6 . 02E+00 

1S.33E+00 

POWER 

( Wh  T \ S ) : 

29.S5E-02 

CHIP  NOS: 

1 6 . S2E+00 

16.  :SE*00 

16. !2E*00 

1 6 . OSE  +00 

16.30E+00 

1 7. 1 1E+00 

POWER 

(WATTS) : 

29.67E-02 

CHIP  N06: 

I6.Ste-00 

I6.50E*00 

16.33E*00 

16.33E-00 

16.S8E+0Q 

18. 13E+00 

POWER 

(WATTS) : 

29.64E-02 

CHIP  NO 7 : 

16.64E+00 

16.42E+00 

16. 12E*00 

16. 12E+0Q 

16.32E+00 

17.31E+00 

POWER 

(WATTS) : 

30.G1E-02 

CHIP  N08: 

16.7  T  E  *00 

16.42E+00 

16.40E*00 

16.40E+00 

16.28E+00 

1 8 . 02E+00 

POWER 

(WATTS) : 

30.12E-02 

CHIP  N09: 

1 6 . 6SE  *00 

16.38E*00 

16.23E*00 

1 6 . 36E+00 

16.41E+00 

1 7.6SE+00 

POWER 

(WATTS) : 

29.37E-02 

HEAT  1 

EXCHANGERS 

TEMPERATURES 

:  RIGHT 

LEFT 

BOTTOM:  ■ 

99. 12E 

-01  99.39E 

-01 

TOP: 

10.33E 

+00  98. OSE 

-01 

BACK  PLANE  TEMPERATURES  ARE 


T<S5>: 
T ( So) : 
T(S7I: 
T( 72) : 
T( 73) : 
T<74>: 
T( 75): 
T c  76) : 
T( 77) : 


1 4 . 67E  *00 
1 S . A ' E  *00 
1S.27E+00 
15.S7E+00 
15.SSE+00 
1 S  .  06E  *00 
! S . SbE  *00 
1S.S6E+00 
IS.4SE+00 


I 


! 


TABLE  2 


TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  0.7  W 


THESE  RESULTS  ARE  STORED  :n  FILE:  D80982100 


AMBIENT  TEMP  WAS:  23  C 

VOLTMETER  READING  WAS:  3.2198  #61 
BATH  TE*P  WAS:  m  C 

THIS  RUN  WAS  RECORDED  WITH  FLOW  VIZ 


ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NOT: 
POWER 

23.3SE-00 
(WATTS) : 

22.5 3E*00 
71 .08E-02 

22.3SE-D0 

22 . 39E*00 

22.24E*00 

24 . 06E+00 

CHIP  N02: 
POWER 

23. 36E*00 
(WATTS) : 

22 . 37E*00 
71.1  IE-02 

22.43E*00 

22.09E+00 

22 . 39E+00 

2S . OSE+OO 

CHIP  N03: 
POWER 

23. I6E»00 
(WATTS) : 

22 . 2G£*00 
70.87E-02 

22.55E*00 

22 . 1 4E+00 

22. 73E*00 

24 .S3E+00 

CHIP  N04: 
POWER 

23.23E*00 
(WATTS) : 

22.42E*Q0 

70.18E-02 

22.31E+00 

21 . 69E+00 

22. 16£*00 

24. 16E+00 

CHIP  NOS: 
POWER 

23. 1 T  E ♦ 0  0 
(WATTS) : 

21  .9SE»00 
70.47E-02 

21 . 94E+00 

21 . 85E*00 

22.37E»00 

24. 19E+00 

CHIP  N06: 
POWER 

22.69E+00 
(WAITS) : 

22. 77E+00 
70.38E-02 

22. 18E*00 

22.21E+00 

22.6QE*O0 

26. 12E*00 

CHIP  NO 7: 
POWER 

23.49E*Q0 

(WATTS): 

22. 76E+00 
71 .26E-02 

22. 1 7E*Q0 

22. 15E+0Q 

22.39E+00 

24.55E*00 

CHIP  N08: 
POWER 

23 . 32E*00 
(WATTS) : 

22.63E*00 
71  .SSE-02 

22.47E+00 

22.S5E+00 

22.47E*00 

26. 1 7E+00 

CHIP  N09: 
POWER 

23.07E-00 
(WATTS) : 

22.54E*0Q 

70.3SE-02 

22.05E+00 

22.38E+00 

22.20E*00 

25.4JE»00 

HEAT  EXCHANGERS  TEMPERATURES: 

RIGHT 

LEFT 

BOTTOM: 

98. 18E-01 

98.46E-01 

TOP: 

98.21E-01 

93.S8E-01 

BACK  PLANE  TEMPERATURES  ARE  : 


T  f  55  : 

17'. 

67E 

♦  00 

T 1 

So) : 

18. 

76E 

♦  00 

T  ( 

57): 

19. 

05E 

♦  00 

T  ( 

72): 

19. 

08E 

♦  00 

T  ( 

73): 

18. 

90E 

♦  00 

T  ( 

74) : 

19. 

1 1E 

♦00 

T  ( 

75): 

19. 

02E 

♦00 

T  ( 

76): 

18. 

79E 

♦00  ' 

T  ( 

77): 

18. 

7 1 E 

♦  00 

TABLE  3 

TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  1.1  W 


THESE  RESULTS  ARE  STORED  IN  TILE:  D8099074S 


AHBIENT  t[;*p  was:  23  C 

VOLTMETER  READING  HAS:  4.1107  *61 
BATH  TEMP  HAS:  10  C 

THIS  RUN  HAS  RECORDED  HITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NOI : 
POHER 

28 . G2E  *00 
(WATTS) : 

27 . 3SE  *00 

1 1  .S7E-01 

27. 1 3E+00 

27. t6E+00 

26.93E»00 

29.81E*00 

CHIP  N02: 
POHER 

28. 73E+00 
(WATTS) : 

27 . 0SE*00 

1 1  .S3E-01 

27. 13£*O0 

26 . 72E*00 

27. 1 1£*00 

31 .43E+00 

CHIP  N03: 
POHER 

28.47E+00 
(WATTS) : 

27. 17£*00 
11  .S4E-01 

27 . 49E+00 

26 . 82E+00 

27.S9E*00 

30.6SE+Q0 

CHIP  NO 4: 
POWER 

28.47E*00 
(WATTS) : 

27. 04£*00 
11  .43E-01 

27 . 1 4E+Q0 

26. 18E+00 

26. 95E*OQ 

30. J3E+00 

CHIP  NOS: 
POWER 

28.21E*00 
(WATTS) : 

26.29£*00 
11 .48E-01 

26.42E<-00 

26 . 30£*00 

27. 03£*00 

30. OOE+OO 

CHIP  NOG: 
POHER 

27.45E+00 
(WATTS) : 

27.90E*00 

11.46E-01 

26 . 98E+00 

27.00E*00 

27. 76E+00 

33.28E+00 

CHIP  N07: 
POHER 

29. 12E*00 
(WATTS) : 

27.77E*00 

11.60E-01 

27. 03E+00 

2S.97E+00 

27.2SE+00 

30 .66E+00 

CHIP  N08: 
POHER 

28. S6E+00 
(WATTS) : 

27.4SE-00 
11 .6SE-01 

27.27E*00 

27 . 40E+00 

27 . 18E*00 

33.20E*00 

CHIP  N09: 
POWER 

28. 13E+00 
(WATTS) : 

27.32E»00 

11.5SE-01 

26.60E+00 

27. 12E+00 

2S.72E-00 

32.07E»00 

$ 

W 

$; 

iV*; 

nil 

js 


M 
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HEAT  EXCHANGERS  TEMPERATURES: 

bottom:  • 

TOP: 


RIGHT  LEFT 

99.44E-01  99.77E-01 

1 0 . 0SE*O0  96.S0E-01 


BACK  PLANE  TEMPERATURES  ARE  : 


T(5S) 
T(S6) 
T  ( S  7  ) 
T  ( 72  > 
T(  73 ) 
T(  74) 
T(7S> 
T  <  76 ) 
T(  77) 


19.?4E*00 
21 .13E+00 
21 . 79E*00 
21 . 42E+00 
21 .  UE*00 
21  . 41 E  *00 
21  .S3E*00 
20 . 93E  «-00 
20 .32E*00 


i 

I 
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TABLE  4 

TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  1.5  W 

THESE  RESULTS  ARE  STORED  IN  FILE:  D8099201S 


AMBIENT  temp  WPS:  23  C 
VOLTMETER  READING  HAS :  *4.5442  *51 

BATH  TEMP  HAS:  'U  C 

THIS  RUN  HAS  RECORDED  WITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NO 1 : 
POWER 

32.9IE+00 
(WATTS) : 

31 . 1 9E*00 
14.77E-01 

3 1 . 03E  »00 

31 ,08E*00 

30 . 78E+Q0 

34.47E*00 

CHIP  N02: 
POWER 

32. 92E*00 
(WATTS) : 

30. 93E+00 
14.77E-01 

31 .02E+00 

30 . 42E*00 

3t .Q7E*00 

36.44E+00 

CHIP  N03: 
POWER 

32. 36E+00 
(WATTS) : 

30.95E*00 

14.72E-01 

31 . J2E*00 

30.5QE+00 

31 .66E+00 

3S.33E*00 

CHIP  NO 4: 
POWER 

32.60E+00 
(WATTS) : 

30.86E*00 

14.S8E-01 

31 . 04E+00 

29 . 86E+00 

30 . 79E+00 

34 .86E*00 

CHIP  N05: 
POWER 

32.28E+00 
(WATTS) : 

29 . 79E+Q0 
14.64E-01 

30 . 1 2E+00 

29.99E+00 

30.81E+00 

34.66E»Q0 

CHIP  nog: 
POWER 

30.86E+00 
(WATTS) : 

31 . 77E+00 
14.62E-01 

30.66E+00 

30 . 72E*00 

31 . 43E*00 

39.92E*00 

CHIP  N07: 
POWER 

33.43E+00 

(WATTS): 

31 .83E*00 
14.8QE-01 

30.89E+0Q 

30.81E*00 

31 . 12E+00 

3S.4SE+00 

CHIP  N08: 
POWER 

32.82E*00 

(WATTS): 

31.45E>00 

14.87E-01 

31 .27E+00 

31 .40E+00 

3! .22E+00 

38. 77E+00 

CHIP  N09: 
POWER 

32.21E»00 

(WATTS): 

31 . 07E+00 
14.74E-01 

30.3GE+00 

30.97E-00 

30.30E+00 

37.2SE*00 

I 
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HEAT  EXCHANGERS  TEMPERATURES: 
bottom: 
top: 


RIGHT  LEFT 

99.49E-01  99.69E-01 

99.69E-01  95.79E-01 


BACK  PLANE  TEMPERATURES  ARE 


T  ( 

SS): 

•> ; 

87E 

♦00 

T  ( 

56): 

23! 

45E 

♦00 

T  <  5  7 ) : 

24. 

33E 

♦00 

T  ( 

72): 

23. 

84E 

♦00 

T  ( 

73): 

23. 

52E 

♦00 

T  ( 

74): 

23. 

83E 

♦00 

T  ( 

7S) : 

24. 

1 OE 

♦00 

T  ( 

76): 

23. 

25E 

♦  00 

T  ( 

77): 

23. 

I OE 

♦00 

S' 

5S 

I 


?• 
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TABLE  5 


TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  1.9  w 


THESE  RESULTS  ARE  STORED  IN  PILE:  D8 1001130 


AMBIENT  TEMP  WAS:  23  C 

VOLTMETER  READING  WAS:  5.2344  #61 
EA Th  TEMP  WAS:  to  C 

THIS  RUN  HAS  RECORDED  WITH  PLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NOl : 
POWER 

37 . 4SE*00 
(WATTS): 

35 . 42E»00 
19. 1 OE -01 

3S.24E*00 

35.29E*00 

34 . 90E*00 

39.63E*0Q 

CHIP  N02: 
POWER 

37 . 4GE+00 
(WATTS): 

34.63E+00 
19. t 1£-0t 

35.00E+00 

34.29E*00 

34 . 99E+00 

42.43E»00 

CHIP  N03: 
POWER 

36.60E+00 
(WATTS) : 

34 . 89£*00 
19.05E-01 

35.33E+00 

34 . 33E*00 

35.51E+00 

40 . 46E*00 

CHIP  N04: 
POWER 

37.22£*00 
(WATTS) : 

34.88E*00 

18.86E-01 

35.28E*00 

33 . 81E+00 

35.0IE*00 

40.20E*00 

CHIP  NOS: 
POWER 

36. 70E*00 
(WATTS) : 

33.55E»00 

18.94E-01 

34.01£*00 

33 . 86E+00 

34.97£«-00 

39. 83E+00 

CHIP  NOG: 
POWER 

35.05E+00 

(WATTS): 

36.S5E+00 

18.90E-01 

34.59£*00 

34.9SE-00 

35.73E+00 

48.38E+00 

CHIP  N07: 
POWER 

3S.33E+00 

(WATTS): 

36. 1 7E*00 
19.1SE-01 

3S. 13E»00 

34.99E+00 

35.34E+00 

40 . 79E*00 

CHIP  N08: 
POWER 

37 . 38E*00 
(WATTS): 

35.8IE+00 

19.23E-01 

3S.42E*00 

35.57E*00 

35.43E+00 

4S.03E+00 

CHIP  N09: 
POWER 

36.47E+00 
(WATTS) : 

3S.22E*00 

19.07E-01 

34.36E+00 

35. 10E*00 

34 . 40E*00 

43. I3E»00 

HEAT  EXCHANGERS  TEMPERATURES:  RIGHT  LEFT 

BOTTOM:  99.S4E-01  10.02E+00 

TOP:  .  95.29E-0!  92.47E-01 


BACK  PLANE  TEMPERATURES  ARE  : 


23.S4E  *00 
25.35E*00 
26.58E+00 
2S.81E+00 
25 . 38E  +00 
2S.72E*00 
26.29E*00 
25.04E+00 
24.8 SE*00 
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TABLE  6 


TEMPERATURE  DATA  FOR  INPUT  POWER  Q^n  =  2.3  W 


THESE  RESULTS  ARE  STORED  IN  FILE:  D81001740 


AMBIENT  TEMP  UAS:  23  C 

VOLTMETER  READING  HAS:  5.8093  *61 
BATH  TEMP  WAS:  10  C 

THIS  RUN  HAS  RECORDED  WITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  N01 : 
POWER 

42. 48E+Q0 
(WATTS) : 

39.94E+00 

23.07E-01 

39. 75E*00 

39.90E*00 

39.39E*00 

45.05E+00 

CHIP  N02: 
POWER 

42.44E+00 
(WATTS) : 

39.2SE*Q0 

23.08E-01 

39.57E+00 

38.83E+00 

39.51E*00 

48.53E*00 

CHIP  N03: 
POWER 

41 .90E+00 
(WATTS) : 

39.S0E+QQ 

23.00E-01 

40.09E*00 

38 . 99E*00 

40.S9E+GO 

46. 1 7E*00 

CHIP  N04: 
POWER 

42.09E*00 

(WATTS): 

39.24E+00 

22.78E-01 

29.82E*00 

38 . 09E*00 

39.48E*00 

45.71E*00 

CHIP  NOS: 
POWER 

4 1 . 49E*00 
(WATTS) : 

37 . 80E*00 
22.88E-01 

38.3SE*00 

38. 1 7E+00 

39.S0E*00 

4S.28E+00 

CHIP  N06: 
POWER 

39.67E*00 

(WATTS): 

41.33E*00 

22.82E-01 

38.98E+00 

39.52E*00 

40 . 48E*00 

56.04E*00 

CHIP  N07: 
POWER 

43.53E+00 

(WATTS): 

40.80E*00 

23.12E-01 

39. 71E*00 

39.58E+00 

39.89E+00 

46.S3E+00 

CHIP  N08: 
POWER 

42.48E>00 

(WATTS): 

40 . 36E*00 
23.22E-01 

40.09E+00 

40.20E+00 

39.57E+00 

51 .52E*00 

CHIP  N09: 
POWER 

41.70E+00 
(WATTS) : 

39.86E*00 

23.03E-01 

38.85E+00 

39.61E*00 

38.95E+00 

49.58E+00 

HEAT  EXCHANGERS  TEMPERATURES:  RIGHT  LEFT 

BOTTOM:  .  1 0 . OOE*00  10.0AE+00 

TOP:  10.1 4E+00  10.Q6E-00 


BACK  PLANE  TEMPERATURES  ARE  : 

T(S5):  25.S9E+00 
T ( 56  J :  27.76E+00 
T(57):  29.42E*00 
T(72):  28.37E*00 
T(731 :  27.87EM10 
T(74>:  28.23E*O0> 

T( 751 :  29.04E*0Q 
T(76):  27. A1E»00 
T <  77 ) :  27 . 1 9E*00  ' 
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TABLE  7 


TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  =  2.7  W 


THESE  RESULTS  ARE  STORED  IN  FILE:  D810II130 


AMBIENT  TEMP  UAS:  23  C 

VOLTMETER  READING  WAS:  5.244  «6t 
BATH  TEMP  WAS:  (0  C 

THIS  RUN  WAS  RECORDED  WITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

BACK 

CHIP  NOt: 
POWER 

46.68E+Q0 

(WATTS): 

43.58E+00 

27.60E-01 

43.37E+00 

43 . G9E+00 

43.06E*Q0 

49.81E-00 

CHIP  N02: 
POWER 

46.2SE*00 
(WATTS) : 

43 . 02£*00 
27.61E-01 

43. 1 4E  *00 

42 . Q4E+00 

42.86E*00 

53.82E*00 

CHIP  N03: 
POWER 

4S. 70E+Q0 
(WATTS) : 

43. 18E*00 
27.S2E-01 

43.85E*00 

42.33E+00 

44.32E*00 

50 . 88E+00 

CHIP  M04: 
POWER 

«6.22E*on 
(WATTS) : 

42.88E-00 

27.26E-01 

43.55E+00 

41 . 56E*00 

43.20E*00 

50.G2E*00 

CHIP  N05: 
POWER 

45.45E*00 
(WATTS) : 

40.31E*00 

27.37E-01 

41 .67E*00 

41.SSE*00 

43. 18E+00 

49.94E+00 

CHIP  N06: 
POWER 

43.55E+00 

(WATTS): 

45. 08E+00 
27.29E-01 

41.98E*00 

43.23E+00 

44.00E+00 

63.01E+00 

CHIP  N07: 
POWER 

47.S6E+00 

(WATTS): 

44.G9E+00 

27.66E-01 

43.43E+00 

43.20E+00 

43.56E+00 

5l.81E»00 

CHIP  N08: 
POWER 

46. 05E+00 
(WATTS) : 

43.92E*00 

27.78E-01 

43.5SE+00 

43. 77E+00 

43.34E-00 

57.24E*00 

CHIP  N09: 
POWER 

44. 78E+00 
(WATTS) : 

43.22E*00 

27.S5E-01 

41 .99E*00 

42.90E+00 

41 .90E*00 

55.35E*00 

HEAT  EXCHANGERS  TEMPERATURES: 
BOTTOM:  • 

TOP: 


BACK  PLANE  TEMPERATURES  ARE 


RIGHT  LEFT 

99-.27E-01  99. 49E-01 

1 0 . QOE+OO  98 . 46E-01 


T  ( 55 1 
T(56) 
T  ( 57 1 
T(  72) 
T(  73 ) 
T(  741 
T(  75 ) 
T(7S) 
T(  77 ) 


26.87E*G0 
29 . 09E*00 
31 .  13E+00 
29.93E+00 
23.33E+0Q 
29. 77E+09* 
30 . 78E*00 
28.82E-00 
28 .57E*00 
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TABLE  8 

TEMPERATURE  DATA  FOR  INPUT  POWER  Qin  -  3.1  W 

THESE  RESULTS  ARE  STORED  IN  TILE:  D8I0HSI5 


AMBIENT  TEMP  WAS:  23  C 

VOLTMETER  READING  UAS:  6.8013  *61 
BATH  TEMP  WAS:  10  C 

THIS  RUN  WAS  RECORDED  WITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 


*  rv*  V./WV  i  m  V.  '  Vi '  v  \r 


CENTER 

TOP 

RIGHT 

LEFT 

bottom 

BACK 

CHIP  NOt: 
POWER 

S0.27£*00 
(WATTS) : 

46.37E*00 
31 .S3E-01 

46.67E*00 

46. 36E+00 

46.34E*00 

54.20E*00 

CHIP  N02: 
POWER 

49. 47E+00 
(WATTS) : 

4S.72E*00 
31 .60E-01 

45. 94E*00 

44.88E*00 

45.65E*00 

58. 34E*00 

CHIP  N03: 
POWER 

48. A9E*00 
(WATTS): 

45 . 75E*00 
31 .S0E-01 

46. 70E*00 

44.97E*00 

47.04E*00 

S4.80E*00 

CHIP  N04: 
POWER 

49. 72E*00 
(WATTS) : 

46.06£*00 

31.19E-01 

46.89E*00 

44.65E+00 

46.51E*00 

55.03E+00 

CHIP  NOS: 
POWER 

48.67E*00 

(WATTS)’. 

43.79E*00 

31.32E-01 

44.88E»00 

44.50E*00 

46. A5E*00 

54. 10E*00 

CHIP  N06: 
POWER 

46.42E*00 

(WATTS): 

48.26E+00 

31.23E-01 

44.90E»00 

46. 1 7E*00 

47.08E*00 

68. 1SE+00 

CHIP  N07: 
POWER 

51 .67E*0G 
(WATTS): 

48.20E-00 

31.66E-01 

46.82E+O0 

46.51E*00 

46.95E*00 

56.86E*00 

CHIP  N08: 
POWER 

49.44E*00 

(WATTS): 

46.99E+00 

31.80E-01 

46.6 1 E  *00 

46.80E*00 

46.32E*00 

62.31E*00 

CHIP  N09: 
POWER 

48.20E+00 
(WATTS) : 

46. !9E*00 
31 .S3E-01 

44.97E*00 

46.00E*00 

45.25E+00 

60.8IE+00 

HEAT  EXCHANGERS  TEMPERATURES: 
bottom: 

TOP: 


BACK  PLANE  TEMPERATURES  ARE  : 


RIGHT  LEFT 

1 0 . 02E  *00  1 0 . 05E»00 

1 0 . 03E*Q0  99.64E-01 


T155) 
T(S6) 
T(S7J 
T(72) 
T(  73 ) 
T(  74  > 
T(  75) 
T ( 76) 
T  <  7  7 ) 


28.S8E«00 
30.93E*00 
33. 10E*00 
31 . 64£*00 
31 .01E*00 
31  .47E*0a 
32.S6E*00 
30 . A6E+00 
30.21E*00 
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TABLE  9 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  0.3  W 


THE  RAW  Emt  DATA  ARE  FROM  THE  FILE: 
THE  POWER  SETTING  PER  CHIP  WAS: 


D8070092S 
0.3  WATTS 


ONET(W) 


Tavg-Ts 


■/.UNC  IN  Nu 


27.38E-02  S6.03E-O1  10.61E+00  20.55E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  104.43E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :179.04E-02 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  M0.76E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.88E-02 

27.38E-02  S7.S0E-01  10.33E+00  20.14E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  107.44E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :l74.32E-02 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  110.98E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.88E-02 

27.29E-02  57 . 70E-Q1  10.28E+00  20.12E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  107.63E-02 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :174.03E~02 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  110.63E-G2 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.88E-02 

27 . 03E-02  55. ISE-Ot  10.S5E+00  20.32E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  102.58E-02 
V.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :182.08E-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  109.23E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.89E-02 

27.14E-02  56.S8E-01  10.42E+00  20.42E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  105.40E-02 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :177.A9E-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  109.86E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.89E-02 

27. 1 1E-02  58.16E-01  10.13E+00  20.00E-0I 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  108.54E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :172.66E-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  109.94E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.89E-02 

27.45E-02  56.19E-01  10.62E+00  20.52E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  104.64E-02 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  : 178. 7 IE-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  111.08E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.88E-02 

27.56E-02  •  58 . OOE-O 1  10.33E+00  20.04E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  108.23E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :173.13E-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  111.75E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.87E-02 

27.33E-02  S7.25E-01  10.37E+00  20.23E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  106.74E-02 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :175.40E-02 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  T  T 0 . 70E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  100.88E-02 
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TABLE  10 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  0.7  W 


THE  RAW  Enf  DATA  ARE  FROM  THE  FILE:  D30701S3S 

THE  POWER  SETTING  PER  CHIP  WAS:  0.7  WATTS 

CHIP  ONET(W)  Tavg-Ts  Nu  7.UNC  IN  Nu 

1  70 . 30E-02  1 1 .51E+00  13.32E+00  13.3SE-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  228.83E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.27E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  304.72E-02 
2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 

2  70.32E-02  11.S1E+00  13.32E+00  1 3 . 35E-0 1 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  228.84E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.24E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  304.80E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 

3  70.09E-02  11.51E+00  13.28E*00  13.3SE-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  228.78E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.44E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  303.80E-02 
2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

4  69.41E-02  1 1 . 36E+00  13.32E+00  13.43E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  225.48E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :883.80E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  300.34E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  t01.07E-02 

5  69.70E-02  11.23E+00  13.S3E+00  13.49E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  222.64E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :893.83E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  301.13E-02 
2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

6  69.6TE-02  11.34E+00  13.38E+00  13.44E-01 

TtMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  224.97E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :885.56E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  301.12E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  1 01  - 07E-02 

7  70.49E-02  11.57E+00  13.28E+00  13.32E-01 

TEMPERATURE -BASED  RAYLEIGH  NUMBER  *  E-6  IS:  230.J9E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :867.69E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  305.75E-02 
2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.05E-02 

8  70.78E-02  •  t1.S2E+00;  13.28E+00  13.30E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  231.24E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :864.18E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  307.15E-02 
2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 

9  70.17E-02  4 1 . 34E+00  13.49E+00  13.43E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  22S.08E-02 

2  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :88S.19E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  303.58E-02 
1  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 
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TABLE  11 

REDUCED  DATA  FOR  INPUT  POWER  Qin  =  1.1  W 


THE  RAW  Eat'  DATA  ARE  FROM  THE  FILE:  D807H035 

THE  POWER  SETTING  PER  CHIP  WAS:  1.1  WATTS 

CHIP  ClNET(W)  Tavg-Ts  Nu  7.UNC  IN  Nu 

1  11.45E-0!  1 6 .67E*00  T  5 . 0 1 E-00  1 1 . 77E-0 1 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  343.22E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :602.83E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  S1S.31E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

2  11. AGE-01  16.59E+00  1S.08E+00  11.78E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:  341.44E-02 

V.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :60S.54E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  51S.06E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

3  11.42E-01  16.56E+00  1S.06E+Q0  11.79E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  340.72E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :606.64E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  5 13. 2 IE-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

A  11.31E-01  16.41E+00  1 5 . 05E+00  11.82E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  337.0SE-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :612.23E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  S07.35E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

5  11.35E-01  16.0SE+00  15.4SE+00  11.89E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  328.29E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :626.l2E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  £-7  IS:  507.3SE-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

6  11.34E-01  15.92E+00  15.55E+00  11.92E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  32S.2SE-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  : 631 . 1 1E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  50S.88E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

7  11.48E-01  16.81E+00  14.93E*Q0  11.74E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  346.68E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :597.78E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  517.48E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.O6E-O2 

8  11.S3E-0I  *  16.76E+00  15.03E+00  11.7SE-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  34S.42E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :599.62E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  519.33E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  !01.0SE-02 

9  11.43E-01  -  16.28E+00  1S.34E+00  11.84E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  333.97E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :617.04E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  512.18E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 
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TABLE  12 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  1.5  W 


THE  RAW  Emt  DATA  ARE  FROM  THE  FILE:  D80720730 

THE  POWER  SETTING  PER  CHIP  WAS:  1.5  WATTS 

CHIP  QNET(W)  Tavo-Ts  Nu  7.UNC  IN  Nu 

1  1 A . 59E-01  19.0SE+0G  16.76E*Q0  11.41E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  401.43E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :S27.67E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  672.71E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.12E-02 

2  14.59E-01  I9.08E+00  16.74E*00  11.40E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  402.2SE-02 

V.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :526.77E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  673.17E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101. T2E-02 

3  14.S4E-01  1 9 . 3SE  *00  16.45E+00  1  1.37E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  409.23E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :519.31E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  673.03E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.1 2E-02 

4  M.40E-01  18. 73E+00  16.83E+00  11.4SE-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  393.3SE-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :536.67E-03 
FLUX  BASED  RAYLEIGH  NUMBER  •  E-7  IS:  661.82E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.13E-02 

5  14.46E-01  18.39E+00  17.20E+00  U.49E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  385.00E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  13  :546.38E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  SS2.05E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.13E-02 

6  14.44E-01  18.59E+00  16.99E+00  1 1 . 47E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  »  E-6  IS:  389.86E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :540.68E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  662.46E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.13E-02 

7  14.62E-01  13.27E+00  16.60E+00  11.38E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  407.23E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  I521.42E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  676.12E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.1  IE-02 

8  14.68E-01  '  1 9 . 48E*00  16.50E+00  11.35E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  412.33E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  : 5 1 6 . 06E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  68O.47E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.1  IE-02 

9  14.56E-01  ’19. 09E+00  16.69E+00  1  1.40E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  402.60E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :S26.39E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  671.79E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.12E-02 
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TABLE  13 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  1.9  W 


THE  RAW  Eat  DATA  ARE  FROM  THE  ^ ILE : 
THE  POWER  SETTING  PER  CHIP  WAS: 


D8072 l 8 i 0 
1.9  WATTS 


QNET(W) 


Tava-Ts 


7.UNC  IN  Nu 


18.99E-01  24.92E+00  16.7SE+00  10.38E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:  56S.77E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :403.27E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  947.69E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

19.00E-0T  24.72E+00  16.89E+00  10.89E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  SG0.03E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :406.S6E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  945.98E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

1 8 . 94E-01  24 . 68E+00  16.86E+00  10.90E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  SS9.04E-Q2 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :407.13E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  942.52E-02 
1  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 

18.75E-01  24.53E+00  16.80E+00  I0.91E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  554.73E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :409.6SE-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  931.79E-02 
7.  UNCERTAINTY  IN  FLUX  8ASED  RAYLEIGH  NUMBER  IS:  101.09E-02 

18.83E-01  23.88E+00  17.32E+00  10.95E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  S36.40E-02 
V.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :420.80E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  928.99 £-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 

18.80E-01  23.83E+00  17.32E+00  10.95E-0! 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  53S.00E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :421.68E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  926.84E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.09E-02 

19.04E-01  25.23E+00  16.S9E+00  10.86E-01 

TEMPERATURE-  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  S74.77E-02 
X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :398.25E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  953.30E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

19. T  2E - 0  T  ‘  2S.20E+00  16.68E+00  10.86E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  573.88E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :398.74E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  957.00E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

18.96E-01  -24 . 47E+00  17.02E+00  10.91E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  SS2.91E-02 
7  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :410.73E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  941.24E-02 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 
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TABLE  14 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  2.3  W 


THE  RAW  Ent'  DATA  ARE  FROM  1  HE  FILE:  D80731030 

THE  POWER  SETTING  PER  CHIP  WAS:  2.3  WATTS 

CHIP  ONET(W)  Tavg-Ts  Nu  7.UNC  IN  Nu 

1  22.85E-01  28.31E+00  17.76E+00  1 0. 71E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  6S8.S3E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :35S.02E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  116. 98E-0 1 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101 . 07E-02 

2  22.86E-01  27 . 97E+00  17.98E+00  10.73E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  648.40E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :359.32E-03 
FLUX  BASED  RAYLEIGH  NUMBER  »-  E-7  IS:  116.60E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

3  22 . 78E-01  27.97E+00  17.93E+00  10.73E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:  648.36E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :3S9.34E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  116. 222-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  1 01 . 07E-02 

4  22.S6E-01  27.37E+00  17.31E+00  10.73E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  645.46E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :360.S9E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  114.97E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 

5  22.65E-01  27 . 1 OE+Ofl  18.38E«-00  10.77E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  622.74E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :370.80E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  114.48E-0I 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 

6  22.60E-01  27 . 1 5E+00  18.31E+00  10.76E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  £-6  IS:  624.26E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  : 370 . 1 0E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  114.28E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.08E-02 

7  22.90E-01  28 . 72E+00  17.55E+00  10.70E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  67I.04E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :349.89E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  1I7.80E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 

8  23.00E-01  ‘  2B.59E+00  17.71E+00  10.70E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  667.0SE-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :351.S1E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  118.13E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  F.AYLEIGH  NUMBER  IS:  101.06E-02 

9  22.80E-01  "  27 . 68E  +  00  18.13E+00  10.74E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:  639.59E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :363.12E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  115. 97E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.07E-02 


TABLE  15 


REDUCED  DATA  FOR  INPUT  POWER  Qin  =  2.7  W 


THE  RAW  Emt  ORTA  ARE  FROM  THE  FILE:  080731766 

THE  POWER  SETTING  PER  CHI?  WPS:  2.7  wATTS 

CHIP  ONET(W)  Tavg-Ts  Nu  7.UNC  IN  Nu 

1  27.37E-01  33 . 28E+00  18.16E+00  10.SSE-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  •*  E-6  IS:  813.30E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :302.08E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  147.70E-01 
X  UNCERTAINTY  IN  FLUX. BASED  RAYLEIGH  NUMBER  IS:  101.05E-02 

2  27.39E-01  32 . 77E+00  18.44E+00  10.S6E-0! 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  796.9SE-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :306.72E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  146.98E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  I01.05E-02 

3  27.30E-01  32 . 53E+00  18.52E+00  10.57E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  789.23E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :308.97E-Q3 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  146.13E-0I 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 

4  27.03E-01  32 . 7SE+00  18.22E+00  10.S6E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  796.28E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :306.91E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  14S.04E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

5  27.1 SE-01  31 . 73E+00  18.87E+00  10.59E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  763.51E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :316.79E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  144.06E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.06E-02 

6  27.07E-01  31.79E+00  18.78E+00  10.S9E-OI 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  765.48E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :316.18E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  143.75E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 

7  27.44E-01  33.82E*00  17.92E+00  10.S3E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-5  IS:  831.14E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :2 97.22E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:  148.91E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.05E-02 

8  27.S6E-01  '  33 . 41E+00  18.21E+00  10.54E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  817.56E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :300.90E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  148.90E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.04E-02 

9  27.32E-01  '32.1 1 E*00  18.77E+00  10.58E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  77S.54E-02 

X  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  : 31 3 . 07E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  145.59E-01 
X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  101.0SE-02 
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1.  Graphical 


resentation  of  Reduced  Data 


The  experimental  surface  temperature  data  were 
plotted  as  a  dimensionless  Nusselt  Number  vs.  both  a 
temperature  based  and  a  flux  based  Rayleigh  Number,  Figures 
4.15  through  4.32. 

Representative  error  bounds  were  plotted  in  Figures 
4.15  and  4.16  for  input  powers  of  1.5  and  1.9  watts.  All 
points,  with  the  exception  of  1.5  watts  were  well  within  the 
error  bounds.  The  point  representing  1.5  watts  was  either 
too  high  on  the  Nusselt  Number  axis  or  too  low  on  the 
Rayleigh  number  axis. 

The  factors  affecting  the  Nusselt  Number  were  h,  L  & 
K.  Of  these  three,  only  h  was  calculated  with  uncertainty 
in  this  experiment.  The  factors  affecting  the  Rayleigh  Num¬ 
bers  were  a,  L,  AT,  Qnet,  v  and  K*  of  these  only  AT  and 
Qnet  had  uncertainties. 

Examination  of  Tables  4  and  12  resulted  in  no  data 
which  appeared  to  be  extreme  in  either  the  raw  or  reduced 
data  sets.  One  possible  explanation  for  this  point  is  the 
changing  flow  patterns  around  this  power  input.  This  may  be 
a  transition  region  from  somewhat  orderly  flow  to  completely 
chaotic  flow.  It  is  recommended  that  further  investigation 
of  this  region  of  power  settings  be  conducted. 

The  equations  for  the  best  fit  lines  are  very  close 
to  each  other  suggesting  that  an  overall  correlation  that 
would  produce  good  approximations  is  possible. 
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TEMP  BASED  RAYLEIGH  NUMBER 

Figure  4.19  Temperature  Based  Rayleigh  Numbers,  Chip  3 

CHIP  NUM8ER  3 


Figure  4.20  Flux  Based  Rayleigh  Numbers,  Chip  3 
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C.  TEMPERATURE  VARIATIONS 


The  temperature  variations  with  time  were  investigated 
and  reported  in  Figures  4.33  through  4.40.  These  results 
were  motivated  by  the  studies  conducted  by  Liu,  Yang  and 
Kelleher  [Liu  et  al.,  1987,  pp.  486-97].  Three  thermocouple 
channels,  numbers  1,  13  and  32  were  scanned  individually  at 
a  rate  of  3  times  per  second  for  a  period  of  200  seconds. 
These  quantitative  results  were  then  plotted.  The  results 
clearly  demonstrate  a  multiple  frequency  content  signal  for 
most  of  the  power  levels  investigated  in  this  experiment. 

An  observation  of  this  experiment  was  that  the  output  of 
thermocouple  #1  remained  steady  with  time  throughout  the 
increase  in  power  levels.  The  temperature  level  for  this 
thermocouple  was  also  generally  higher  than  for  the  other 
two  thermocouples  sampled. 

A  possible  explanation  for  these  results  is  that  the 
cold,  stagnant  fluid  layer  created  by  the  bottom  heat 
exchanger  prevented  the  penetration  of  warmer  fluid  below 
the  bottom  row  of  components.  This  is  also  supported  by 
flow  visualization. 

One  exception  is  at  3.1  watts.  As  seen  in  Figure  4.40, 
there  is  a  large  fluctuation  of  temperature  with  time.  A 
possible  explanation  for  this  fluctuation  may  be  that  the 
increased  energy  acquired  by  the  fluid  at  this  power  setting 
enabled  the  flow  to  force  itself  into  the  stagnant  region  at 
the  bottom  of  the  chamber. 
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V.  RECOMMENDATIONS 


There  is  a  great  deal  more  experimentation  to  be 


accomplished  with  this  component  configuration. 


Three 


possible  variations  to  the  chamber  are  : 


-  Placing  a  spacer  in  the  chamber  to  reduce  the  chamber 
width.  This  would  act  in  the  same  manner  as  a  shrouding 
wall  and  would  affect  the  heat  transfer  characteristics. 


-  Setting  the  bottom  heat  exchanger  to  a  temperature  other 
than  10°C.  The  energy  expended  to  maintain  the  bottom 
heat  exchanger  at  10°C  is  considerable.  If  reasonable 
components  temperatures  can  be  realized  by  maintaining 
the  bottom  heat  exchanger  at  a  higher  temperature  level 
then  the  energy  consumption  of  the  system  would  be 
reduced . 


Powering  different  configurations  other  than  the  entire 
array.  This  could  possibly  offer  some  explanation  of  the 
flow  patterns  observed  in  the  chamber. 


Other  possible  improvements  to  software  and  hardware 


-  Producing  a  plotting  program  that  would  read  the  data 
directly  from  the  data  acquisition  unit  for  plots  of 
Nusselt  vs.  Rayleigh  number. 


-  Implementation  of  measurement  devices  that  have  greater 
sensitivity  to  minor  fluctuations  in  temperature 
variation. 


-  Development  of  a  program  that  will  perform  frequency 
analysis  on  the  data  recorded  by  the  program  Fastscan. 
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APPENDIX  A 


Throughout  this  appendix,  the  sample  calculations  will 


be  based  on  component  number  one  with  a  power  input  of  0.71 


watts. 


1.  Determination  of  Input  Power 


Figure  A. 1  Electrical  Network  of  Power  Input 


From  Equation  (3.4): 


Qin  =  ((3.218  -  0.5278)  (0.5278) )/2 


IV.' • 


7<V» 


2.  Estimation  of  Conduction  Loss 


Equation  3.3  requires  Ts,  Tj-,  and  Rcon.  Ts  and  Tb  were 
read  directly  from  thermocouples  placed  immediately  behind 
the  heaters  and  on  the  back  of  the  board,  respectively. 
Rcon  was  calculated  by  using  the  constriction  resistance 
formulas  developed  by  Mikic  and  repeated  in  [Kraus  and  Bar- 
Cohen,  1983,  pp.  76-9].  The  conduction  losses  for  components 
3 ,  5  and  7  were  then  averaged  to  yield  a  total  conduction 

loss. 

3 .  Calculation  of  Nusselt  and  Rayleigh  Numbers 

The  calculations  of  the  convection  coefficient  from 
Equation  (3.5)  was: 


Qnet  “  ^(Tavg  “  Tc) 


net 

h  "  A  (T  -T  ) 
avg  c 


.703 


(5.76  x  10-4) (11.51) 


=  106.04  W/m2  °K 


From  Equation  (3.10): 


T  +  T 
avg  c 


Tfilm  “ 


Tfilm  =  15 . 75°C 


Using  the  relationships  of  [3-M  Manual,  1985,  pp.  9-60],  the 
physical  properties  of  FC-75  were: 

kf  =  (.65  -  7.8947  X  10"4  *  Tfilm)(0.1) 

=  0.0637  W/m°C 

p  =  1000(1.825  -  0.00246  *  Tfilm) 

=  1823.21  Kg/m3 

Cp  =  4180(. 241111  +  3.7037  x  10~4  *  Tfilm) 

=  1032.23  J/Kg°C 


v  =  10~6(1.4074  - 

2. 

964  X  lO-2  *  Tfilm 

+  3.8018 

X 

10-4  *  Tfilm 

-  2.7308 

X 

10"6  *  Tfiim 

+  8.1679 

X 

10”9  *  Tfiim) 

v  =  1.024  x  10“6  m2/sec 

_  _ .00246 _ 

B  “  1.825  -  .00246  *  Tfilm 

8  =  1.377  X  10“3 
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Then  from  Equation  (3.1): 


106 .04  ( .008) 
.0637 


=  13.32 


.0637 

1823.2  *  1032.2 


=  3.38  x  10-8  m2/sec 


1.024  x  10 
3.38  x  10-8 


=  30.3 


Grt  = 


gBL  (T  -  T  ) 
y  avg  c 


9.81(1.377  x  10~3)  (.008)3  (11.51) 
(1.024  x  10-6)2 


=  7.592  X  10‘ 


i 


From  Equation  (3.11),  the  temperature  based  Rayleigh  number 
was: 

Ra^  =  Gr^  *  Pr 
=  2.3  x  106 

From  Equation  (3.12),  the  flux  based  Rayleigh  number  was: 
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'SQloss  =  °*055  *  Qioss 


■  .000159  W 


Therefore, 


SQnet  "  -0071 


The  uncertainty  for  the  other  variables  used  were: 


6Tsink  =  5T 


=  0.1  K 


4Tfilm  ~  ( ^Tsink  +  ®avg) 1/2 


=  (.012  +  o.l2)1/2 


.1414  K 


TDelt  *  (sr|vg  +  «T?sink)  i/2 

=  (O.l2  +  .oi2)1/2 


=  .1414  K 


* 
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wnet. 2  ,  xPelt'2 

Qnet  TDelt 
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APPENDIX  C 


SOFTWARE 


1.  Acquire 


10 

20 

30 

40 


41 

42 
SO 
60 
70 
80 
90 
100 
110 
120 
130 
140 
ISO 
160 
170 
180 
190 
200 
210 
211 
220 
230 
240 
2S0 
260 
261 
270 


280 

290 

291 
300 
310 
320 
330 
340 
3S0 
360 
370 
380 
390 
400 

1  08" 


410 

420 

430 

440 

4S0 

460 

470 

480 

490 


S00 

S10 

S20 

S30 


!  FILE  ACQUIRE 
I 

!  EDITED  BT  LI  T  J  BENEDICT.  USN 

j 

1  READ  FILE  "READ_ME" 
i 

COM  /Co/  D( 7) 

DIM  Emf  t  76 1  .Power ( 9 ) , T ( 76 1 
! 

1C0RRELATI0N  FACTORS  TO  CONVERT  EMF  TO  DEGREES  CELCIUS 
DATA  0.10086091 . 2S727 . 9. - 7G734S . 8 . 7802SS96 , 

DATA  -924  7486S89.6.98E1 1 . -2 . 66E  1  3 . 3 . 94E  1'4 
! 

READ  D<«> 

Rp-2.000 
PRINTER  IS  701 
BEEP 


INPUT  "ENTER  THE  INPUT  MODE:  O-SYS.  1-FILE". In 
I 

IF  Im-1  THEN 
BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  FILE  TO  BE  READ" .Oldf i leS 

i 

PRINT  USING  " 1SX , * " THESE  RESULTS  ARE  STORED  IN  FILE  : 

i 

ELSE 

BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  NEH  FILE" .Newf i leS 
PRINT  USING  " 1  OX. " " THESE  RESULTS  ARE  STORED  IN  FILE: 
END  IF 


PRINT 

PRINT 

i 

INPUT  "ENTER  THE  AMBIANT  TEMP" . TemoS 

PRINT  USING  "1SX, ’"AMBIENT  TEMP  HAS:  "" . 10A" : TenpS 


"".lOA'lOldfileS 


*",10A":Neuf ileS 


INPUT  "FLOU  V 117  Y /N" , AnsS 
INPUT  "VOLTMETER  READING". VS 

PRINT  USING  "1SX,""V0LTMETER  READING  HAS:  "".IOA’jVS 

t 

INPUT  "ENTER  THE  BATH  TEMP’.BS 

PRINT  USING  "1SX,""BATH  TEMP  HAS  :"",10A";BS 

IF  AnsS-"Y"  THEN  PRINT  USING  "1SX.*"THIS  RUN  HAS  RECORDED  HITH  FLOH  VIZ"". 


IF  In- 1  THEN  ASSIGN  «File  TO  OldfileS 
PRINT 

i 

IF  In-0  THEN 
CREATE  BDAT  Newfl leS.S 
ASSIGN  «Flle  10  NewfileS 
END  IF 

i 

!  READ  DATA 

i 

IF  In-0  THEN 

OUTPUT  709; "AR  AFOO  AL79* 


87 


hh 


600 

OUTPUT  722: *F1  R1  T1  ZO 

FLO’ 

sso 

i 

S60 

FOR  1*0  TO  76 

S70 

OUTPUT  709:*AS* 

S80 

WAIT  1 

SSO 

ENTER  722 :E»T ( I ) 

600 

BEEP 

610 

NEXT  I 

620 

1 

610 

OUTPUT  AFi 1« : Emf  («) 

600 

1 

6S0 

ELSE 

660 

ENTER  $Fiie;E«f<«) 

670 

END  IF 

680 

i 

690 

OUTPUT  709; *TD* 

700 

i 

710 

FOR  1-0  TO  60 

720 

Su»*0 

730 

FOR  J-1  TO  7 

700 

Su«*Sun*D( J)*Emf ( I ) *  J 

7S0 

NEXT  J 

760 

T ( I ) *Su« 

770 

NEXT  I 

780 

i 

790 

FOR  1*71  TO  76 

800 

Su«*0 

810 

FOR  J*1  TO  7 

820 

Sun*Sum+D( J)#Enf(I)‘J 

830 

NEXT  J 

800 

T ( I ) -Sun 

3S0 

NEXT  I 

860 

1 

870 

1  POWER  CALCULATIONS 

880 

! 

890 

j- 1 

900 

Vol  t-E»>T  (61  ) 

910 

( 

920 

FOR  1-62  TO  70 

930 

Power (J) -Emf (I)*(Volt-Ewf(I 

900 

J-J*1 

9S0 

NEXT  I 

960 

1 

970 

BEEP 

980 

BEEP 

990 

PRINT  USING  *  MX,  "ALL 

TEMP 

1000 

i 

1010 

PRINT 

1020 

i 

1030 

PRINT  USING  *12X. "CENTER 

1030  PRINT  USING  * I2X, "CENTER  TOP  RIGHT  LEFT  BOTTOM 

BACK'" 

1000  PRINT 

I0S0  PRINT  USING  *  IX,  "CHIP  N01 :  "  ,S< ZZ. DDE . 2X >  .2X.ZZ. DDE* : 1(0) . T (I ) . T (2) , T1 
. T(0> ,T<S> 

1060  PRINT  USING  *5X,  "POWER  (WATTS):  "  , ZZ . DDE* :Power ( 1 ) 

1061  PRINT 

1070  PRINT  USING  *1X."CHIP  N02:  ** ,5( ZZ . DDE ,2X ) , 2X. ZZ . DDE* ; T (6) . T ( 7 ) . T ( 8* , T ( 
,1(10). Kill 

1080  PRINT  USING  *SX, "POWER  (WATTS):  " , ZZ . DDE* -.Power ( 2) 

1081  PRINT 
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CalcDiel 

•  FILE  NAME  :  CalcDiel 

i 

!  EDITED  BY  :  Lt.  T  BENEDICT.  USN 

i 

!  THIS  PROGRAN  ANALYSES  THE  DATA  READ  FROM 
!  A  DATA  FILE  DESINATED  BY  THE  OPERATOR.  IT 
!  REDUCES  THE  DATA  TO  CALCUALTIONS  OF  NET 
!  POWER,  RAYLEIGH  NUMBER  AND  NUSELT  NUMBER. 

! 

CUM  /Co/  D(/> 

! 

DIM  Emf(70t, Power ( 9).T(60).Tava<9).Ts(9) 

DIM  Tfilm(9).Qnet(9),H(9).K(9).Rho(9).Cp(9> 

DIM  N<9) ,Nu(9) ,Ra(9> .Delt<9> .Alfa (9) ,Pr<9) 

DIM  Gr(9) ,Beta(9) ,Dpou(9) ,Dts(9) 

j 

!  CORRLEATION  FACTORS  TO  CONVERT  Emf  TO  REGREES  CELCIUS. 

DATA  0.10096091 ,25727 . 9 , -76734S . 8 . 7802SS96 . 

DATA  -9247486589, 6. 98E 1 1 , -2.66E13.3.94E1 4 

i 

READ  DC*) 

i  , 

!  RESISTANCE  OF  THE  PRECISION  RESISTOR  IS:  17. 

Rp-2.00 

i 

PRINTER  IS  701 

BEEP 

BEEP 

i 

INPUT  “ENTER  THE  NAME  OF  THE  FILE  CONTAINING  DATA" ,01df ileS 

PRINT  USING  “1SX,**THE  RAW  Emf  DATA  ARE  FROM  THE  FILE:  ** , 10A* :01df i ieS 

INPUT  "ENTER  THE  POWER  SETTING  “, Power* 

PRINT  USING  *  1 7X,  “  THE  POWER  SETTING  PER  CHIP  HAS:  ** , 1 0A“ ;PouerS 

PRINT 

BEEP 

BEEP 

ASSIGN  AF  i  le  TO  OldfileS 
ENTER  8File,‘Emf <*) 

i 

i  CONVERT  Emf  TO  DEGREES  CELCIUS 

FOR  1-0  TO  60 

Sum-0 

FOR  J-0  TO  7 
Sun-Sum+D ( J ) *Emf ( I ) " J 
NEXT  J 
T ( I ) -Sum 
NEXT  I 

i 

i  CONVERT  Emf  TO  POWER 
J-1 

Vol t-Emf <61 ) 

FOR  1-62  TO  70 

Power! J) -Emf ( I )*(Voit-Emf ( I) ) /Rp 
Dpow ( J ) - . 0 1 *Po«e  r  <  J ) 

J-J+J 
NEXT  I 

i 

!  AREA  OF  THE  BLOCK  FACES 
Acen-1 .92E-4 


640  Ar lg* 1 .44E-4 
650  Atop*4.8E-5 

660  Abot-4.8E-S 
670  Atot-S . 76E-4 
660  ! 

690  ! CALCULATE  THE  AVERAGE  TEMPERATURES  OF  THE  BLOCK  FACES 

700  ! 

710  Tavg( t ) - ( T ( 0 ) *Acen+T ( 1 >*Atop  +  T<2)*Ang+T(3)»Alef  +  T(4)*Abot )/Atot 

720  Tavg(2)*(T<6)*Acen+T(7)*Atop  +  T(3)*-Ario+T(9l*Aiet+T(10)*Abot)/'Atot 

730  Tavg(3>*( T ( 12)*Ac<?n+T( 13>*Atop+T( 14)*Arig+T( IS)*Alef *T< 1 6 ) *Abot ) /A  tot 

7 40  Tavgl 4  > ■ ( T ( 1 8 ) *Acen*  r(19)»Atop+T!20)*Arig*T(21>*Alef+T(22)*Abjt)/Atot 
7S0  Tavg(5) *<  T(24 )*Acen+T (2S )*Atop+T <26)*Arig+T(27)*Alef+T <28 )*Abot ) /Atot 
760  Tavg(6)-(T<30)*Acen*T(31 )  *A top  +  T ( 32 i *Ar ig+  T ( 33 ) *A] ef  +  T ( 34 )  *-Abot ) /A tot 

770  Tavg<7)-< T ( 36 > *Acen+T ( 37 ) «Atop ♦ T ( 38 > *Ar xg+T ( 39 ) *Alef +  T ( 40 > *Abo t ) /A  to t 

780  Tavg (8)*(T(42>*Acen+T(43) * A  top ♦T(44)»Aria+T(4S) *A leT+T(46)*Abot)/Atot 

790  Tavg(9>-(Tl48)*Acen+T(49)*Atop+T(50)*Arig+T(51 )*Alef»T(S2)*Abot)/Atot 

800  ! 

810  !  UNCERTAINTY  OF  THE  THERMOCOUPLE  READINGS 

820  Dt* . 1 
830  Dtavg*.1 
840  ! 

3S0  !  RESISTANCE  OF  PLEXIGLASS  AND  UNCERTAINTY 

860  Rc-764.6S 

880  ! 

890  !  CHIP  BACK  SURFACE  TEMPERATURES 

900  Ts( 1 >*T(5) 

910  Ts(2)-T<m 

920  Ts (3) *T ( 1 7 ) 

930  Ts<  4 ) »T( 23 ) 

940  Ts (S ) *T( 29) 

9S0  Ts(S)-T(3S) 

960  Ts( 7 ) *T( 41 ) 

970  Ts(8)-T ( 47 ) 

980  Ts( 9) *T( S3) 

990  Tssum-0 
1000  FOR  J-1  TO  9 
1010  Tssum*Tssum+Ts< J) 

1020  NEXT  J 
1030  ! 

1040  Tsavg*Tssum/9 
10S0  ! 

1060  !  CONDUCTION  LOSS  CALCULATION  AND  UNCERTAINTY 

1070  Qloss3-(T(17)-T(S6))/Rc 

1080  Qloss5* ( T( 29) -T(55) ) /Rc 
1090  0 loss7- ( T ( 4 1 )-T(S4) )/Rc 

1  100  Qloss- ( Qloss3+Q.loss5+Qloss7 ) /3 
1110  Dqloss*(QIoss)*.022 
1120  ! 

1130  !  AVERAGE  SINK  TEMPERATURE  CALCULATION 

1140  Tsinlr-(T(57)+T(58)+T(59)  +  T<60)  )/4 

1150  Dtsmif  *(Dt/T(58) ) 

1160  L-8.0E-3 

1170  ! 

1180  PRINT  USING  ’IX,  "CHIP  ONET(H)  Tavg-Ts  Nu 

7.UNC  IN  Nu  ‘  ,  1  0A" 

1190  PRINT 
1200  ! 

1210  ! 

1220  !  CALCULATION  OF  NET  POWER.  Nu.  Ra.  AND  UNCERTAINTIES 


I 

I 

I 

r» 


m 


M 

fell 


M 
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m 


1230  ' 

1240  FOR  J-1  TO  9 
124!  PRINTER  IS  1 
1242  PRINT  “J” 

1244  PRINT  J 

1245  PRINTER  IS  701 

1250  ! 

1260  !  CALULATIQN  OF  Onet 

1270  One t ( J) -Power ( J i -Q loss 
1280  Dqnet'(DpowlJ) "2+Dqioss‘21 ‘ .5 
1290  ! 

1300  !  CALCULATION  OF  Tfilm 

1310  Tf i lm( J ) • ( Tavg( J ) +  Ts ink > /2 
1320  Dtf lira* (Dts ink ‘2+Dtavq‘21 * .5 
1330  ! 

1340  !  CALCULATION  OF  A  DELTA  TEMPERATURE 

1350  Delt( J)*Tavg( J)-Tsink 
1360  Ddelt*(Dtavg‘2+Dtsink ‘2) ‘ .5 
1370  ! 

1380  !  CALCULATION  OF  CONVECTION  COEFFICIENT 

1390  H( J)-Onet( J)/(Atot*Delt( J)> 

1400  Dh'H(J)*( (Dqnet/Qnetl J)) *2+iDdelt/Delt( J) ) *  2 ) * . 5 

1410  ! 

1420  !  CALCULATION  OF  CONDUCTION  COEFFICIENT 

1430  K<  J ) -  C  .65-7.89474E-4*Tf  i  imU)  )/10 
1440  ! 

1450  !  CALCULATION  OF  DENSITY 

1460  Rho(J)-(1 .825-. 00246*Tf1lm(J)) *1000 

1470  ! 

1480  !  CALCULATION  OF  SPECIFIC  HEAT 

1490  Cp( J) -( .241 1 1 1*3. 7U37E-4*Tf llm(J) 1*4180 

1500  ! 

1510  !  CALCULATION  OF  VISCOSITY 

1520  N ( J ) —  1  . 4074-2. 964E-2*Tf i Iffl< J) *3.8018E-4*Tf ilm( J) *2-2.7308E-6*Tf ilm(J) ‘3+8. 
l679E-9*Tf llml J) '4 
1530  N( J ) *N( J )* 1 . E-6 
1540  ! 

1550  !  CALCULATION  OF  BETA 

1560  Beta(J)-. 00246/(1 .825-. A3246*Tfilm(J)) 

1570  ! 

1580  !  CALCULATION  OF  ALPHA 

1590  Alfa(J>*K(J)/( Rho ( J ) *Cp ( J ) ) 

1600  ! 

1610  !  CALCULATION  OF  PRANDTL  NUM8ER  AND  UNCERTAINTY 

1620  Pr( J)'N( J)/Alfa( J) 

1630  ! 

1640  !  CALCULATION  OF  NUSSELT  NUMBER 

1650  Nu(J)-H(J)*L/KCJ) 

1660  Dnu'Nu ( J) * ( Dh/H( J) ) 

1670  ! 

1680  !  CALCULATION  OF  PERCENT  UNCERTAINTY  IN  NUSSELT  NUMBER 

1690  Pernu* ( Dnu/Nu< J) ) * 1 00 
1700  ! 

1710  !  CALCULATION  OF  GRASHOF  NUMBER  AND  UNCERTAINTY 

1720  Gr( J) *9.81*Beta( J)*(L‘3)*Del t ( J) /N( J) "2 
1730  Dgr-Gr(J)*(Ddelt/Deit(jn 
1740  ! 

1750  !  CALCULATION  OF  RAYLEIGH  NUMBER  AND  UNCERTAINTY 

1760  Ra(J)-Gr(J)*Pr( J)*1 .E-6 
1770  Dra-Ra<J>*<Dgr/Gr(J>)  ■ 

1  780  ! 


1794  !  CALCULATION  OF  FLUX  BASED  RAYLEIGH  NUMBER 

1800  ! 

1310  Raf(J)*' (9.81«Beta(  J)*L'4*0neUJ) ) / (K f J ) *N( J) *A1 f a (J I *Atot )) »  I  .E-7 
1320  Draf  *Raf  t  J )  *  (  Dqne  t/Cinet  l  J  )  I 
1830  ! 

1340  1  PERCENT  UNCERTAINTY  IN  RAYLEIGH  NUMBERS 

I860  Per ra* ( Dra/Ra l J ) ) *  1 00 
1860  Pe r raf * ( Dra t /Rat ( J ) ) *  1 00 
1870  ! 

1880  PRINT  USING  " 1 X . D . 1 X . 4( 7X . ZZ . ODE , 1  * ; J . Onet ( J ) . Del t ( J ) , Nu ( J ) .Pernu 
1830  ! 

1910  PRINT  USING  “6X, ""TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:  " " .3D.DDE" : 
Rat  J) 

1920  PRINT  USING  "6X.**7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  I 
S  :  " * . ODD . DDE" : Per  ra 

1930  PRTNT  USING  “6X.’"FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:  ’ " . 3D . DDE’ : Raf ( J) 

1940  PRINT  USING  "6X,""7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:  “ " .  ODD. 

DDE" :Perraf 

1950  PRINT 

1960  NEXT  J 

1970  ASSIGN  SFile  TO  * 

1980  END 


m 


I 

|| 


w 


i  j,'  t.*  <.»  ir:.- *  Jt*. 


3 .  FastScan 
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\  PROGRAM  FASTSCAN 

!  EDITED  Bt  LT.  T  J  BENEDICT 
! 

•  THIS  PROGRAM  SCANS  3  CHANNELS  FOR  VOLTAGE  MEASUREMENTS 

Ipass-593 

Pass-0 

N-0 

DIM  T1 (S99> .VI  (?) .Yt (V39) 

DIM  T2<  S99 ) . V2( 2 ) . Y2<  S93 ) 

DIM  T3 (S99) . V3<2) , Y31S99) 

CLEAR  709 
CLEAR  722 

i 

BEEP 

PRINTER  IS  701 
PEEP 

INPUT  ’ENTER  THE  FIRST  FILE  NAME:  '.Newfilett 
INPUT  ’ENTER  THE  SECOND  FILE  NAME:  *.Newfile2S 
INPUT  ’ENTER  THE  THIRD  FILE  NAME:  *.Neufile3* 

INPUT  ’ENTER  THE  VOLTMETER  READING:  *.V« 

PRINT  USING  ’ 1SX , ’ *  THESE  RESUL’S  ARE  NOW  STORED  ON  DISK  FASTSCAN' 
PRINT 

PRINT  USING  *2SX.”FILE:  ”  . 1 OA’ ;Ne«f i let S 
PRINT 

PRINT  USING  *2SX. ’ ’FILE :  ”  . 1 0A* ;Newf  i  le2S 
PRINT 

PRINT  USING  *2SX. ’’FILE :  * ’ , 1 0A* ;Ne«f i le3* 

PRINT 
WAIT  1 
BEEP 

OUTPUT  709: *AE I  * 

WAIT  2 
BEEP 

OUTPUT  722; "TA  FI  R1  PO  ZO  1STI  SOT  1STN* 

! 

!  LOOP  NUMBER  ONE 

I 

OUTPUT  709; *AF0  ALO’ 

OUTPUT  709;’AS’ 

BEEP 

T  imedatel-TIMEDATE 
FOR  Jj-0  TO  Ipass 
OUTPUT  722: *T3* 

ENTER  722:V1<»> 

T  1  (Pass  >  -  VI ( ( ) 

Pass-Pass*  1 
N-N*  I 
NEXT  J) 

T  i«edate2*TIMEDATE 

Tot a  1 1 met * T i»edate2-T imedatel 

OUTPUT  722;*AC0’ 

Pas*-0 

I 

1  LOOP  NUMBER  TWO 

i 

OUTPUT  709; ’API 2  AL12* 

OUTPUT  709:’AS* 

BEEP 


”,  10A- 


■"-.’VW 


FOR  J  i - 0  TO  Ipass 
OUTPUT  722 :  *  T3" 

ENTER  722:V2<*> 

T2<Pass>  -V2(  1  ) 

Pass-Pass*  I 
NEXT  J, 

OUTPUT  722; "AC12* 

Pass-0 

i 

!  LOOP  NUMBER  THREE 

i 

OUTPUT  709;  'AF3T  AL31* 

OUTPUT  709; *  AS* 

BEEP 

BEEP 

BEEP 

FOR  J  i  *0  TO  Ipass 
OUTPUT  722; ' T3* 

ENTER  722:V3< •) 

T3(Passi-V3< I ) 

Pass-Pass*  t 
NEXT  J; 

!  END  LOOPS 

PRINT  USING  ‘1SX.*‘THE  TOTAL  TIME  ELAPSED  WAS  (SECONDS ):** ,2X. ( DDD. DDE > * 

1 1 i me  1 

PRINT 

PRINT  USING  '  f  5X, *  *  THE  TOTAL  NUMBER  OF  SCANS  WAS  :  *  *  .  2X,  ( DDDD.  D.2X )  *  ,*N 

PRINT 

PRINT  USING  *1SX.**THE  VOLTMETER  READING  WAS:  *M0A.*;V* 

PRINTER  IS  1 
BEEP 

i 

j  TRANSFER  FIRST  SCAN  DATA 
! 

CREATE  BDAT  NeufiletS.20 
ASSIGN  «F,le  TO  NewfilelS 
OUTPUT  £F i lei T 1 <  * ) 

FOR  I i-O  TO  Ipass 

THI.  )-.1008609!*2S727.9«T1(Ij)-76739S.8*T1(Ij),2*78002S56«T1(Ii)'3 
NEXT  Ii 
! 

!  TRANSFER  SECOND  SCAN  DATA 
! 

CREATE  BDAT  NeUf,le2S.20 
ASSIGN  «F i )e  TO  NewT  i  le2S 
OUTPUT  <*F  j  le;  T2<  *  » 

FOR  I i -0  TO  Ipass 

T2< 1 1 ) * . 1 008609 T *25727 ,9«T2( Ii)-76734S.8*T2(Ii)'2* 78002SS6* T2(  1 1 ) ‘ 3 
•EXT  Ii 
I 

«  TRANSFER  THIRD  SCON  DATA 
! 

CREATE  BDAT  Newfile3S,20 
ASSIGN  $File  TO  New{ile3S 
OUTPUT  «Fi le;T3(«) 

FOR  I i -0  TO  Ipass 

T3( 1 1 )* . I00B6091 *2S727.9#T3( Ii)-76734E.8»T3(Ii) ‘ 2* 78002556*  T3( 1 1 ) *  3 
NEXT  I. 


4.  Plot 


FILE  NAME:  PLOT 
EDITED  BY  Lt  T  J  BENEDICT 

THIS  PROGRAM  PLOTS  THE  DATA  ACQUIRED  BY 
PROGRAM  -FASTSCAN* . 

80  PRINTER  IS  70S 

90  BEEP 

100  Xmin*0 
1 1  0  Xmax*200 
120  BEEP 

130  INPUT  “ENTER  MINIMUM  AND  MAXIMUM  Y-VALUES* .Ynwn.Ymax 

140  BEEP 

ISO  Xstep*20 

160  BEEP 

170  Ystep”. 2 

180  BEEP 

190  PRINT  “IN:SP1 : IP  2000.2000.8000,7000:* 

200  PRINT  “SC  0.100.0, 100; TL  2,01* 

210  SFx”100/(Xmax-Xmn) 

220  Sfy*100/<Ymax-Ymin) 

230  PRINT  “PU  0.0  PD" 

240  FOR  Xa“Xnin  TO  Xmax  STEP  Xstep 

250  X-(Xa-Xnun)*Sfx 

260  PRINT  “PA“;X,*,0;  XT;* 

270  NEXT  Xa 

280  PRINT  “PA  100.0;PU:“ 

290  PRINT  *PU  PA  0.0  PD* 

300  FOR  Ya-Ymn  TO  Y«ax  STEP  Ystep 

310  Y«( Ya-Ymin)*Sfy 

320  PRINT  *PA  0.*;Y,*YT* 

330  NEXT  Ya 

340  PRINT  “PA  0.100  TL  0  2* 

350  FOR  Xa-Xnin  TO  X«iax  STEP  Xstep 

360  X-(Xa-Xnun)*Sfx 

370  PRINT  “PA“;X,*,100;  XT* 

330  NEXT  Xa 

390  PRINT  “PA  100,100  PU  PA  100.0  PD* 

400  FOR  Ya*Ymin  TO  Ynax  STEP  Ystep 

410  Y*( Ya-Ymin)*Sty 

420  PRINT  *PD  PA  100,*,Y,*YT* 

430  NEXT  Ya 

440  PRINT  “PA  100.100  PU“ 

4S0  PRINT  *PA  0.-2  SR  1 .5.2* 

460  FOR  Xa*Xmxn  TO  Xmax  STEP  Xstep 
470  X*  ( Xa-Xn»  in )  *SFx 

480  PRINT  *PA“;X.,*.0;* 

490  PRINT  “CP  -2,-1 ;LB*;Xa;** 

500  NEXT  Xa 

510  PRINT  *PU  PA  0.0* 

520  FOR  Ya-Ynun  TO  Ymax  STEP  Ystep 
530  IF  ABStYaXI  .E-5  THEN  Ya-0 
540  Y-(Ya-Ymin)*Sfy 
550  PRINT  *PA  0.*;Y,** 

560  PRINT  *CP  -S.-.2S;LB*;Ya,*** 

570  NEXT  Ya 
580  BEEP 

590  IdI-0 

600  IF  Idl-0  THEN 


JWPJ^  U'WV. _ 


.  -  -  AV.y  .v.v  ww^. 


610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1021 

1030 

1040 

1050 


‘Time  ( sec  ) 1 


“CP  o.o- 


"ENTER  THE  NAME  OF  THE  DATA  FILE", D_f i leS 
TO  D  files 


BEEP 
X  labels* 

BEEP 

Ylabel$““Temperature  (C)“ 

PRINT  "SR  1.S.2:PU  PA  50. -10  CP*:-LEN1XlabelS)/2:*0:LB*:XlabelS: 
PRINT  "PA  -11.50  CP  0,“:-LEN(YlabeIS)/2*5/6;“DI  0.1  IlB” ; Y laoe 1$ ; 
END  IF 
PRINT 
BEEP 
INPUT 

ASSIGN  SFile 
BEEP 
(1d*0 
BEEP 

Npa i rs-600 
BEEP 

PRINTER  IS  1 
Sym*  1 

PRINTER  IS  70S 
PRINT  ”PU  Dr" 

IF  Sym- 1  THEN  PRINT 
IF  Sym-2  THEN  PRINT 
IF  Sym-3  THEN  PRINT 
IF  Nd> 1  THEN 
FOR  1-1  TO  (Md-1) 

ENTER  9File;Xa.Ya 
NEXT  I 
END  IF 

FOR  Xa-0  TO  199  STEP  .3333333 
ENTER  @F  i  lei Ya 

Ya- .  1 008609 1 +25727. 9*Ya-767345. 8*Ya‘2+78002556*Ya‘3 
X- ( Xa-Xmin ) *Sf  x 
Y-( Ya-Ymin)«Sfy 
IF  Sym>3  THEN  PRINT  "SM" 

IF  Sym<4  THEN  PRINT  "SR  1.4. 2. 4* 

PRINT  *PA* .X.Y.’PD* 

IF  Sym>3  THEN  PRINT  *SR  1.2. 1.6' 

IF  Sym-4  THEN  PRINT  "UC2. 4. 99.0. -8. -4. 0.0. 8.4.0. ;* 

IF  Sym-5  THEN  PRINT  "UC3. 0,99. -3. -6. -3, 6. 3. 6, 3. -6:' 

IF  Sym-6  THEN  PRINT  "UCO .5 .3 .99 . 3. -8 . -6 .0 .3 .8: ' 

IF  Sym-7  THEN  PRINT  'UCO. -5. 3. 99,  -3. 8. 6. 0,-3. -8;' 

NEXT  Xa 
PRINT  "PU" 

BEEP 

ASSIGN  <?FiIe  TO  * 

END 


*SM." 

*Sli+* 
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